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ABSTRACT 


The behavior of capillary jets of non-Newtonian 
inelastic liquids and of viscoelastic fluids was invest- 
igated experimentally. Data were obtained through high 
speed still photographs and motion pictures taken along 
the length of the liquid column. The" program concentrated 
on profile decay measurements, which refers to the change 
m jet diameter with axial distance due to stress relaxation, 
and on jet breakup studies. 

The visccinelastic liquids all possessed yield stresses 
and were comprised of a wide variety of organic and in- 
organic materials. The less gellified fluids formed capillary 
jets which propagated disturbances m the same manner as 
for a Newtonian jet, that is, as an exponentially growing 
wave of constant wave length. Even the most highly gellified 
fluids could be readily atomized. The behavior of this 
class of non-Newtonian liquids could be explained if it was, 
assumed that the time required to reform the gel structure 
outside the capillary was large compared to the jet breakup 
time. Thus, the 'apparent viscosity possessed at the nozzle 
exit is substantially maintained throughout the lifetime 
of the jet. Using t}iis model, the wavelength and breakup 
length data for the viscoinelastic jets were m fair 
agreement with that predicted for a Newtonian jet of the 
same apparent viscosity. 


xvi 



The viscoelastic test liquids were composed of aqueous 
solutions of 'drag reducing' polymers. These formed capillary 
jets which differed considerably m appearance from Newtonian 
and non-Newtonian inelastic jets. Disturbances were propag- 
ated as a non-linear series of droplets connected by threads 
which thinned with distance and ultimately broke. Even m 
the most dilute systems studied , the initial linear wave 
propagation was arrested and gradually evolved into the 
droplet-thread configuration which characterizes these 
fluids. Within the linear range of wave propagation, the 
data for dilute solutions confirms the prediction of 
stability theory that a viscoelastic jet without structure 
is less stable than a Newtonian jet of the same zero shear 
viscosity. 

High speed motion pictures taken along the length 
of viscoelastic jets show that large scale axial deformations 
occur continually. This shows the jet to be under tension 
since these strains require the existence of axial normal 
forces for their continuance. 


xvii 
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CHAPTER 1 
INTRODUCTION 


1.1 PRELIMINARY REMARKS 

The atomization of a liquid is a phenomena of consider- 
able practical importance and examples of it may be found 
in such diverse fields as spray drying, fuel combustion 
and liquid-liquid extraction. It usually occurs under 
turbulent flow cohditions and is one of the most difficult 
problems m fluid dynamics to treat analytically. Test 
data are correlated empirically. For this reason, the low 
speed breakup of a capillary jet of liquid is often studied 
as an aid towards understanding atomization. Newtonian jet 
breakup is one of the few flow geometries where linear 
stability theory has been able to accurately predict 
observed results. 

Newtonian jets have been studied extensively, both 
experimentally and theoretically. "In recent years, the 
use of non-Newtonian additives has substantially increased 
and yet there is very little published information available 
on the behavior of non-Newtonian jets. This research 
represents an introduction to this important area. 

The term 'non-Newtonian liquid' encompasses fluids 
with a wide range of properties. Though their classification 
into groups of similar characteristics is somewhat arbitrary, 
two categories may be broadly identified- 1. those materials 
which, when m the melt state or m solution, are found to 
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exhibit both viscous and elastic properties, commonly known 
as viscoelastic fluids, and 2. those materials which are 
non-Newtonian but do not possess measureable elastic 
effects. These subdivisions will be used m discussing the 
experimental results for capillary jets of non-Newtonian 
liquids „ 

A natural question arises as to whether a viscoelastic 
jet is more or less stable as compared to a Newtonian jet. 

In other geometries, the question remains unresolved as 
data are available to support both conclusions. A 
theoretical analysis by Middleman (58) using a linearized, 
three constant Oldroyd model (66) shows a viscoelastic jet 
to have a higher growth rate of axisymmetric disturbances 
than a Newtonian fluid of the same zero shear viscosity. 

If the growth rate is assumed inversely proportional to 
the breakup length, an assumption which is satisfactory 
with Newtonian liquids, then it can be concluded that a 
viscoelastic jet will exhibit a shorter breakup length. 
Kroesser and Middleman (39) confirmed this result for 
capillary jets of polyisobutylene in tetralm. 

1.2 AIM AND SCOPE OF THESIS 

Because of the lack of published experimental data 
on the behavior of non-Newtonian jets, there is a need to 
broadly define the breakup characteristics of capillary 
jets produced by various classes of non-Newtonian fluids. 
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A complete description of the experimental apparatus 
and procedures used in this research is given m chapter 3. 
Each chapter presenting the test results also includes an 
abbreviated version of the techniques used. 

Some pertinent physical properties and shear stress- 

shear rate measurements of the test fluids are presented 

m chapter 4 . The available rheological, instruments also 
♦ 

determined yield stresses (the minimum shear stress 
required for flow) on the vxscomelastic liquids but could 
not obtain normal stresses with the viscoelastic fluids. 

Chapter 5 considers the rate of relaxation of stresse 
m an unconfmed liquid column due to the abrupt removal 
of the capillary wall. This phenomena occurs m the region 
immediately outside of the nozzle. The problem is solved 
for the case of a dilute viscoelastic fluid and the predict 
ions are compared to experimental data taken with 
Newtonian, non-Newtonian inelastic and viscoelastic liquids 

The sixth chapter investigates the breakup of a 
laminar capillary jet of a non-Newtonian fluid. High speed 
still and motion picture photographs vividly illustrate 
the different characteristics of Newtonian, viscomelastic 
and viscoelastic jets. A theoretical analysis due to my 
colleague. Dr. J. Yerushalmx (29,94) is included for 
comparision with experimental data. 
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Chapter 7 closes this work with a summary of the 
results and the conclusions drawn from them. 

i 

Experimental data, which cannot be obtained from the 
tables and figures m the main text of the thesis, are 
included m the appendix. The computer,. programs used for 
data reduction are straightforward and have been omitted. 

1.3 THE NOTION OF HYDRODYNAMIC STABILITY 

A fluid flowing in a given geometry is subject to a 
variety of natural disturbances arising from sources such 
as burrs on the walls of the container or external 
vibrations from the environment. No matter how carefully 
the flow system is isolated, there will always be disorders 
on a microscopic level. The question arises as to whether 
these small perturbations will dampen with time or will 
amplify m such a way as to create a new disturbed flow 
pattern. Such patterns are known to exist, the most 
familiar being the transition from laminar to turbulent 
flow. If the disturbance dampens with time, as m laminar 
flow, the system is stable; if it amplifies, the flow is 
unstable . 

Since a given set of flow conditions defines a 

stable regime, while another set can produce an unstable 

% 

state, some criteria must exist that provide a condition 
of neutral or marginal stability. The theory of 
hydrodynamic stability is concerned with the location of 
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this locus. Results are often stated m terms of a dimen- 
sionless parameter whose magnitude must be below some 
critical value to assure stability. Typical 'examples of 
disturbed or secondary flows are given m Table 1-1. 


Table 1-1 . Examples of Secondary Flows 


Mode of Flow Pattern of Instability 

v 


Natural 

Convection 


Hexagonal Benard Cells 


Cylindrical Taylor Vortices 
Couette Flow 

Pipeline Flew Random Eddies 


Critical Parameter 

Rayleigh or 
Marangom Number 

Taylor Number 
Reynolds' Number 


The breakup of a capillary jet is always unstable and, 
m this instance, stability is spoken of on a relative basis 
between different fluids. 

An elementary example will be chosen to illustrate the 
small perturbation technique used m analyzing the effects 
of an infinitesimal disturbance superimposed upon the main 
flow. In the usual summation notation with Cartesian 
coordinates, the equations of continuity and motion can be 
written as 




ax* 


- 0 


■at 


y- 

* *** 


_L 2lP_ 

e axL 


-u v v. 


(i.D 

( 1 . 2 ) 


where U is the kinematic viscosity and V is the Laplacian. 
Body forces have been assumed to be negligible. 
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At: steady state, in the absence of oscillations, the 
velocity and pressure are functions of position only. 

Vl - vL(*Oj P« p(x K ) ( 1 . 3 ) 

For small disturbances superimposed upon the main flow, the 

* 

velocity and pressure may be expressed as 


Vl- 

Vl 4- 

e'V , t0 4 

6 * V™ + ... . 

(1.4) 

p = 

P + 

e'p 0 ' + 

iP* .... 

(1.5) 


where £ is a small parameter and 


Vi«V t (x*,V) 5 P«P(X K >t) 


( 1 . 6 ) 


When equations (1.4) and (1.5) are substituted in the 
equations of motion and continuity, a partial differential 
equation results which can be grouped into similar powers 
of £ . The series of equations is expressed as 



.clYs + W; d-Yv.- - 


St 



(1.7) 


SXv 


at 


T7 SVL + 

V i ^6 


v 0 ) AYs. = 


JO 


J_ dE 
e dx t 


2 (l) 

v 7 V, 


(li8) 


etc. 

If infinitesimally small disturbances are assumed, it is 
usual to neglect terms of £ beyond £ and also to omit 
product terms which are second order in the disturbance. 
With these simplifications, the series (1.8) reduces to 

r 



(1.9) 



( 


This equation is often put into dimensionless form m 
order to more clearly display the controlling parameters. 

By choosing a form of the stream function - ^ it may be 
possible to reduce the partial differential equation to an 
ordinary differential equation. Equation (1.9) admits 
solutions which are exponential in rime s the dependence 
on the space variables- } and Xj comes from the 
physical knowledge of the problem. Thus the disturbance 
can be written as 

oct 

'V = PCx^x^x^e d.io) 

The governing equations together with the boundary 
conditions define a characteristic value problem with c< 
as parameter. If all the characteristic values of o< have 
negative real parts, the disturbances will dampen with 
time and the flow is stable. If just one of the character- 
istic values is positive and real, the flow is unstable. 

Extensive treatments of the linearized theory of 
hydrodynamic stability may be found in the monographs of 
Lm (42) and Chandresakar (8) . 

The disturbances cannot continue to grow exponent- 
ially without bound and the problem must be treated on the 
basis of non-linear stability theory once the disturbances 
are finite. Early contributions have come from Meksyn 
and Stuart (47), and Stuart (82). More recently, 

Eckhaus (18) has summarized his research m a monograph. 
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1.4 PROPERTIES OP NON-NEWTONIAN FLUIDS 

A detailed discussion of the classification and 
properties of non-Newtonian fluids is presented and is 
intended to aid the reader who may not be familiar with 
this area of fluid dynamics. In examining the behavior 
of non-Newtonian fluids, it is useful to group -together 
fluids of similar rheological proper tae s as a means of 
achieving some sort of order. The classification given 

S 

m table 1-2 was adopted from that given by Metzner (51) . 

'Examples of typical flow curves exhibited by each of the 

fluid categories are shown m figure 1 - 1 . 

Table 1-2 » Classification of fluid behavior. 

A. Purely viscous fluids- the shear rate depends only upon 
imposed shear stress. 


where TT , ^ and A are the shear 
stress, the apparent viscosity and 
the rate of deformation respectively. 

1. Newtonians = constant 

2. Non-Newtonian : v? = (I, II, III) or 4 = A ( X ) where 

I, II, III are: the principal invariants of the rate 

of deformation tensor (22). This category includes 
pseudoplastic and Bingham plastic fluids m which ^ 
is a decreasing function of the invariants and also 
embodies dilatant fluids in which **? is an 
increasing function of the invariants . 

B. Time dependent inelastic fluids- the shear rate depends 
upon both the imposed shear stress and the duration of 
the stress* 

X = 7 A J 7 = 7 (I, II, III, t) 

This subdivision describes materials which! are part of 
the groups known as thixotropic and rheopectic fluids. 
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Figure 1-1. Fluid flow curves 

(arithmetic coordinates) 

Pseudoplastic fluids exhibit a decreasing slope as 
the shear rate is increased which is synonymous with a 
decreasing apparent viscosity while dilatant liquids show 
the opposite effect. Bingham plastics do not flow until 
the yield stress is exceeded and then display the linear 
relation between shear stress and shear rate that is common 
to Newtonian liquids. Viscometer measurements on a 
thixotropic fluid would reveal hysteresis loops of the type 
Shown in figure 1-1. Curve 1 illustrates the response of 
a fluid whose yield stress is unaffected' by the shear 
history while curve 2 depicts a liquid m which the recovery 
is slower and the yield value is not constant. This effect 
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xs reversxble xn that suffxcient rest brings about a 
recovery of structure. Rheopectxc fluxds would also 
exhibxt hysteresis loops of simxlar form but wxth the 
arrows drawn opposxte to that shown. 

The following subsectxons summarize the more « 
commonly accepted mxcroscopxc explanations for the 
macroscopic behavior of non-Newtonxan liquids. The 
constitutxve equations given relate to simple shear in 
whxeh there is only one non-zero velocxty component, 
which varies in a direction perpendicular to the velocity 
vector. 

1.4-1 BINGHAM PLASTICS 

A Bingham plastic xs assumed to possess a three 
dxmensional structure whxeh has sufficxent rigidxty to 
resist shearing stresses that are less than the yxeld 
stress. Once the yield value is exceeded, the structure 
breaks down and flows as a Newtonian liquid. It obeys 
the following relationshxps % 

- 'tj =9^, A V2 . ^ 

7 = r l2 .>ry 

where tT is the yield stress and A -j_ 2 the constant 
shear rate for a fluxd subjected to the shearing stress- 


(1.13) 

(1.13a) 

(1.13b) 
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While no real fluid behaves m exactly the manner 
described/ many materials do approximate Bingham plastic 
behavior. Some of the more common examples are drilling 
mud/ chalk suspensions and sewer sludge. 

1.4-2 PSEUDOPLASTIC FLUIDS 

This important category describes the behavior of 
many commonly used non-Newtonian fluids. Pseudoplasti- 
city occurs m solutions oi : asymmetric particle's which are 
loosely bound to one another by secondary forces to form 
a randomly entangled network. Under shear, the particles 
tend to orient -jthemselves so that their major axep are 
parallel to the direction of shear. This behavior would 
also be displayed by systems of solvated particles whose 
degree of solvation was a function of the rate of shear. 
The apparent viscosity decreases with increasing rate of 
shear. At very low shear rates, there is little alignment 
and the structure exhibits Newtonian behavior with a zero 
shear viscosity- . When the shear rate is very high, 

the particles are fully extended and Newtonian behavior is 
-again evidenced with an infinite shear viscosity- 

The change m particle configuration with shearing 
stress cannot occur instantly and pseudoplastic fluids 
must be regarded as those m which the time for structural 
rearrangement is small with regard to the time constant of 
the experiment. 

The following three constitutive equations for 
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pseudoplastic flow represent the more popular ones of the 
many purely viscous relations proposed (3) . None of them 
are able to describe stress relaxation or normal stress 
phenomena, two types of behavior common to viscoelastic 

V 

fluids. 

1. Power law 

a. / "t| s K (1*14) 

K/A l ,r (1.14a) 

b. tw= P _ j/KVoV' for tubular flow (1.15) 

ZJL K ~o^ J 

where XL/ is the wall shear stress m a tube of diaraeter-2a 0 
and length- Q for a fluid flowing with an average velocity- 
V Q under an imposed pressure gradient- AP. 

The power law is a two constant empirical equation 
which many fluids obey over moderate ranges of shear rates . 

Xt fails to describe observed behavior at low shear rates 
in that it predicts an infinite viscosity at zero shear rate. 

2. Ellis model 


A e = <■$•>"] 


(1.16) 

(1.16a) 


This three constant empirical model reduces to 
Newtonian behavior with constant viscosity- 7^ 0 at low shear 
rates and approaches power law behavior at high shear rates. 

3. Eyring-Powell model 

-1 

smh B & ^2 


7 = + (70- 


(1.17) 
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The Eyrmg-Powell fluid, a three constant semi- 
theoretical model, reduces to Newtonian behavior at both 
low and high shear rates while predicting pseudoplastic 
behavior at intermediate shear rates. Its use has been 
limited because it results m forms that are difficult 
to solve analytically o 

1.4- 3 DILATANT FLUIDS 

Dilatant fluids display an apparent viscosity which 
increases with increasing shear rate. This behavaor usually 
occurs in concentrated slurries in which there is msuffic-' 
ient liquid present to enable the particles to flow freely 
past one another. Metzner and Whitlock (56) have compiled 
a number of systems which exhibit dilatancy. 

1.4- 4 TIME DEPENDENT NON-NEWTONIAN FLUIDS 

In materials exhibiting time dependent behavior, the 
shear stress (at a given shear rate) changes with the 
duration of the shear. Those fluids in which the shear stress 
increases with time are called rheopectic and those in 
which a decrease occurs are called thixotropic. 

The mechanism responsible for thixotropy is possibly 
similar to pseudoplasticity. In any unsteady flow, from a 
Lagrangian viewpoint, the orientation of asymmetrical 
molecules and particles or the breakdown of solvated masses 
cannot occur instantaneously with respect to time. 
Pseudoplastic behavior may' be considered as a form of 
thixotropy which has a time constant much smaller than the 
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time scale of the experiment. Carver and Van Wazer (7) 
have shown that certain fluids may require times of several 
minutes or hours before the structure is completely regained. 
Several authors, including Fischer, Bauer and Wiberly (21) 
present evidence that this formation time is also proport- 
ional to the level of shear to which the material was 
previously subjected. Thixotropy occurs m many practical 
materials such as paints, food products and some polymeric 
solutions . 

Rheopexy may be considered as being caused by the 
same phenomena as diiatancy but with time dependent effects 
large enough to be detected. It is a relatively rare 
occurrence but has been observed m certain sols and 
bentonite clay suspensions. 

1.4-5 VISCOELASTIC FLUIDS 

High molecular weight polymers dissolved m solution 
or alone as a melt form a class of non-=Newtonian behavior - 
known as viscoelasticity. These materials differ from 
elastic solids m that they flow under an imposed pressure 
gradient and, unlike purely viscous fluids, they display 
elastic properties such as recoil and stress relaxation. 

Such responses are related to the ability of the polymer 
chain to assume different spatial configurations under 
shear flow. A description of some characteristic 
viscoelastic phenomena follows. 
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A. RECOIL 

If the shearing stresses imposed upon a viscoelastic 
fluid under flow are rapidly removed, the fluid will undergo 
a partial recovery of strain. This may be observed exper- 
imentally by following the movement of tracer particles m 
the solution (40) 6r the deflection of an unconstrained 

t 

rotor m a concentric cylinder apparatus (7). 

The viscoelastic fluid may be said to possess a 
'partial memory' of its rest state. An elastic solid has a 
'perfect memory' and returns completely to its rest state 
when the deforming stress is removed. 

B. JET SWELLING 

An unconfined liquid jet emerging from a cylindrical 
tube can expand to a diameter considerably larger than the 
nozzle diameter, as shown m figure 1-2. 



Moderate shear rates High shear rates 


Figure 1-2. Jet swelling for a viscoelastic fluid.. 

The magnitude of the effect is dependent upon the particular 
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Newtonian liquid Viscoelastic liquid 

Figure 1-4. Pressure distribution m a parallel plate 
shear field. 


D P STRESS RELAXATION 

If a Newtonian liquid is subjected to a flow history 
and then held at constant deformation, the shearing 
stresses will be observed to relax instantaneously, even 
though the strain persists. In a perfectly elastic solid, 
the stress caused by the deformation does not relax at 
all. Between these two extremes are viscoelastic fluids 
which undergo stress relaxation at a finite rate. 

Stress relaxation experiments have been used to 
evaluate the constants and predictions of proposed 
constitutive equations (20) . 

E. DRAG REDUCTION 

The addition of certain polymeric additives m 
concentrations as low as 10-50 parts per million can reduce 
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turbulent pressure drop by as much as 50% (23,36,76,84). 

Seyer and Metzner (78) report that for solutions exhibiting 
drag reduction, the flow may be transitional up to Reynolds 
numbers of about 105. 

F c ELONGAT IONAL VISCOSITY 

Viscoelastic theories predict that, under tension, 
a viscoelastic fluid will exhibit an elongational viscosity 
which increases with shear rate (10,43,92). This is m 
contrast - to the shear viscosity which is a decreasing 
function of shear rate. Experimental data are sparse but 
Ballman (1) reported a constant elongational viscosity for 
polystyrene. The ability of polymers to be drawn into 
long fibers has been related to this predicted property 
(44,64,95) . 

1.4-6 CONSTITUTIVE RELATIONS FOR VISCOELASTIC FLUIDS 

Ideally, the constitutive equation for a viscoelastic 
fluid should be complete enough to describe all of the 
previously described phenomena yet be of such a form as to 
yield tractable solutions from the equations of motion. 

It must obey certain invariance principles (22) since 
natural laws are independent of any coordinate transformation. 
In addition, the response of the material to a given set 
of causes must be indifferent to the coordinate system 
chosen. Both requirements may^be satisfied by casting the 
constitutive equation m a tensorial form and formulating 
the relation m a coordinate system that moves and deforms 
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with the continuum (65) . 

Some general statements can be made concerning what 
the constitutive relation for a viscoelastic fluid should 
contain if it is to predict observed phenomena. 

a) Experimental data show that m the limit of zero 
velocity gradient, the shear dependent viscosity tends 
towards a limiting value. The viscous part of the complex 
viscosity, determined from transient response experiments 
also approaches the same limiting value. Therefore, the 
constitutive equation should also contain a zero shear 
viscosity- T^o - 

b) The constitutive equation should contain a 
characteristic time(s) which can be identified with stress 
relaxation phenomena. 

c) The constitutive equation should correlate steady 
shear as well as oscillatory data. 

d) Normal stresses m steady shear should be predict- 
able from the model. 

e) The shear viscosity, shearing stress and normal 

stress should be functions of the shear rate and, in simple 

shear, be of the experimentally observed form shown in 

figure 1-5. The slope £ 3i should be unity at low? 

d log A 12 

shear rat^ * and less than one at intermediate to high shear 

rates. Similarly, the slope d lo 9 ( T n~ X22^ is initially 

<3 log A 12 

two and is reduced below that value at higher shear rates. 
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Figure 1-5. Experimentally observed response of 
’ viscoelastic fluids. 

Constitutive equations for viscoelastic fluids may 
be classified as being either of the 'rate' type or of the 
' integral * type . 

Rate equations are characterized by time derivatives 
of the stress and deformation* Linear relatiohs may be 
written m the form 


Pt - 

Q> 



(1.18) 

where 

P - 

1 + + 

A &L 

^ Dt 2 " ^ ^ A w 

(1.18a) 

and 

<R- 

7- ( 

pe + /'"ft* ^ 

(1.18b) 


P and Q are linear, time differential operators which 
satisfy the principle of material objectivity, while Aj 
and JU-i are relaxation and retardation times, respectively. 

2? ( ) is a convected time derivative such as the Oldroyd 
(65) or Jaumann (67) derivative. The published experimental 
evidence does not favor a particular derivative over another. 
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( 1 . 21 ) 

( 1 . 22 ) 


where 




Equation (1.22), being a linear formulation, possesses the 
same limitations as do linear rate equations. While it is 
able to predict normal stresses and stress relaxation 
phenomena, the viscosity is constant. If the relaxation 
function is taken to be dependent upon the rate of deform- 
ation tensor, a shear dependent viscosity will be ■'obtained . 
Several non- linear integral theories are summarized m 
Middleman (59) . 

Metzner, White and Denn (54,55) have considered 
criteria for selecting viscoelastic constitutive equations 
which are applicable to solving practical engineering 
problems. The Deborah nomber, defined as the ratio of the 
fluid's relaxation time to the deformation time, is used 
as a natural ordering parameter to determine the appro- 
priate type of relation. 

Several review articles are available ($,45,81) which 
compare the predictive ability of rate and integral 
constitutive equations with experimental data. The compa- 
rison does not argue for e Lther rate or integral constitu- 
tive equations. However, Spriggs, Huppler and Bird conclude 
than selected constitutive equations from each type show 
considerable promise as rheological equations of state for 
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use xn fluxd dynamical calculatxons . 

While the preceding theories have treated the 
material as a continuum, molecular models have been success- 
fully developed which predict non-Newtonian behavior as a 
function of shear rate m terms of known molecular 
parameters. A comprehensive review of several molecular 
theories and comparisons with experimental data may be 
found m Middleman (5§) . 

1.5 A DESCRIPTION OF NEWTONIAN JET BREAKUP 

When liquid flows within a capillary, the confining 

r 

walls serve to retard the fluid movement and thereby give 
rise to shearing stresses. For Newtonian liquids m 
fully developed flow, this results m the well known 
parabolic velocity profile. Upon emerging from the nozzle 
into ambient air, the liquid column is no longer confined 
and undergoes a stress relaxation at a rate determined by 
the flow velocity and the rheological properties of the 
fluid. At low velocities, the surroundings exert no 
influence on the liquid column. The appearance of an 
unconfmed, laminar Newtonian jet is represented m figure 
1-6. It is assumed m this representation that the fluid 
velocity is sufficient to render gravity effects neglig- 
lble and thajt the time required for) relaxation is short 
compared to the lifetime of the jet. 
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Nozzle diameter = 2a 0 Breakup length-L 

Relaxed diameter = 2a Wave length- 6 

Figure 1-6. Breakup of a laminar capillary jet 
of a Newtonian liquid. 


Superimposed upon the liquid column are small disturbances 
from sources such as microscopic burrs within the capillary 
or external vibrations transmitted through the nozzle. 

Surface tension forces alone act to return the locally 
disturbed surface to its equilibrium position but overshoot 
the mark. A larger restorative force results from the over- 
shoot which then produces an even greater deviation from 
equilibrium. The disturbance thus propagates as an 
axially symmetric wave of increasing amplitude and fixed 
wavelength until surface tension forces pinch off the 
column to form individual droplets. The intact portion of 
the liquid column is known as the breakup length. 

The breakup length of a capillary jet of a Newtonian 
liquid as a function of velocity has the form shown m figure 


1-7. 
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and the breakup length again increases. This should not be 
mterperted that turbulence enhances stability since, as 
shown m figure 1-8, the breakup time past point D is a 
decreasing function of jet velocity. 


Breakup 

time 



Figure 1-8. Breakup time of a Newtonian capillary jet 
as a function of jet velocity 
(points correspond to figure 1-7) <> 
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CHAPTER 2 
STATE OF THE ART 


This chapter will review the published literature 
dealing with the breakup of a low speed, capillary jet of 
liquid moving through ambient air. In order to keep this 
survey to a reasonable size, the related stability topics 
of a liquid jet in liquid surroundings and the breakup of 

i 

a turbulent jet have been omitted. When background 
material is needed for completeness, efforts were made to 
cite a recent review article rather than list the numerous 
original contributors. 


2.1 STRESS RELAXATION 

A fluid which emerges from a capillary into 
ambient air suffers an abrupt removal of the confining 
wall. This results m a relaxation of all prior stresses. 
If the relaxation is complete within the lifetime of the 
jet, the velocity profile rearranges from the original, 
radially varying shape to a flat profile. Continuity 
arguements show that these changes must also cause the jet 


diamete 1-2300 to differ from the nozzle diameter-2a 0 and 
00 

the ratio is termed JX= CWa.o • 

For capillary jeis of Newtonian liquids, it has been 


found experimentally that the diameter ratio is a function 
of the capillary Reynolds number-Re 0 alone for fully 
developed flow (26,32,60,62). Gavis and Modan (26) report 
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that a nozzle length to diameter ratio- i(/2a 0 greater than 
32 is needed to produce consistent results for % . At 
Reynolds numbers below 16, Newtonian jets expand due to 
viscous forces (%> I ) while at Re Q <16 contraction occurs 
because of inertial forces. Gavis (25) summarized the data 
of Middleman and Gavis (60) as 


%= 1.10 


0.23e 


-60/(Re o ) 


3/2 


( 2 . 1 ) 


Goren and Wronski (32) did not find any upper limit to 
the diameter ratio at the lower Reynolds number range. 

Only at Reynolds numbers greater than 300 was the lower 
theoretical limit of 'X- - f^/z approached (35). 

Goren (30) considered the rate at which the stresses 
relax in a Newtonian jet. He used a boundary layer 
approach, valid for small axial distances where the 
periphery cjf the liquid column does not interact with the 
core, and showed that the jet diameter varies as the cube 
root of the axial distance. 


(a Q -a) /a Q 


0.703 (z/a 0 Re 0 ) 


1/3 


( 2 . 2 ) 


At large axial distances, where the boundary layer 
assumptions are invalid, Goren and Wronski employ a 
linear perturbation technique and show that the relaxation 
follows an exponential decay of the form 

— hjjz/aoo 

n e 



a 


i 


(2.3) 
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The decay rate-h is a function of the jet Reynolds number 
and a surface tension parameter. 

They experimentally verified the predicted dependence 
for small and large axial distances. For capillary Reynolds 
numbers below 100, the jet diameter reached 99% of its 
final value within 1.5 jet diameters. Middleman (57) 
obtained similar data and found that, at Reynolds numbers 
up to 1000, relaxation is 99% complete within 3 jet 
diameters. Thus, for Newtonian liquids, the relaxation 
length may be neglected m comparison to the breakup 
length over the usual range of Reynolds numbers encountered. 

The relaxation of stresses outside the capillary for 
a non-Newtonian jet involves considerations of the effect 
of a shear dependent viscosity and, m the case of 
viscoelastic liquids, the presence of normal stresses. 

Powell (69) and Middleman and Gavis (61,62) investigated 
stress relaxation in capillary jets of viscoelastic 
liquids. They found that a nozzle length to diameter 
ratio of &/Z.C i 0 ■= 40-50 produced consistent results for ~)C , 
which is similar to the Newtonian value. However, Metzner 
et. al,(52) measured a small but definite effect on ~X, when 
d/za 0 was increased to 126 from 97. These conflicting 
conclusions regarding the entrance length required for 
fully developed flow remain to be resolved. The data from 
all investigators show that, at low shear rates an expansion 
is obtained, X > 1, but as the shear rate is increased con- 
traction always occurs. 
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No experiments have been reported on the direct 
measurement of the rate of stress relaxation m non- 
Newtonian jets. In the experiments just reviewed for the 
effect of shear rate on the diameter ratio, the authors 
make no comment on the axial distance required to achieve 
constant values for % . This may be construed as meaning 

that the axial distance was of the same order as for 
Newtonian liquids. However, the data of other investi- 
gators with capillary jets of viscoelastic liquids 
indicate that the relaxation length may be considerable. 
Gill and Gavis (27) imposed artifical disturbances upon 
the motion of carboxymethyl cellulose jets and measured 
the nodal lengths of the resulting transverse wave. 

Starting from a model of a wave propagating on a string 
of slowly varying tension they obtained, after consider- 
able linearization, an expression which permitted the 
calculation of the fluid's relaxation time from their data. 
Values of 0.32 seconds and 0.27 seconds were calculated 
for solutions of 1.4% SCMC and 1.7% SCMC, respectively. 

This means that a slowly decaying axial tension exists over 
the lifetime of the jet. Goren and Gavis (31) have 
developed a theoretical analysis for this geometry. By 
assuming an exponential decay of the axial normal stress, 
a solution was obtained for the nodal points of the wave 
as a function of the axial distance. Middleman (59) 
summarizes experimental results obtained by Middleman and 
Rautenbach for artificially disturbed SCMC jets. Relaxation 
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times on the order of 10 seconds were calculated and, 
while this value differs by an order of magnitude from those 
reported by Gill and Gavis, it also shows that the relax- 

i 

ation length j is considerably longer than for Newtonian jets. 
Additional evidence m this direction comes from Gaskins 
and Phillipoff (-24) who measured the ratio of the jet 
velocity to the average' capillary velocity-V^/Vo and found 
it to be less for viscoelastic jets than for Newtonian jets. 
They attribute this to the presence of accumulated elastic 
energy which absorbs some of the kinetic energy. Further 
experiments appear to be necessary to clarify the phenomena 
of stress relaxation with capillary jets of non-Newtonian 
liquids . 

Mass and momentum balances have been used (22,25,26, 

60,62,79) to relate the change m jet diameter from that 

of the capillary to the forces causing the change, i.e. 

the velocity profile rearrangement, the normal stress 

relaxation and surface tension forces. Control surfaces 

are taken at the capillary exit and at a point far enough 

downstream for the jet to be completely relaxed except for 

an everpresent axial surface tension force. A properly 

formulated balance must reproduce the experimentally observed 

facts that, for Newtonian jets, the diameter ratio is a 

strong function of the Reynolds number and the Weber number 

has little, if any, effect. Contraction results for Re 0 <* 16 ; 

1 

expansion occurs for Re 0 > 16 and, at Reynolds numbers ^greater 
than 300, the diameter ratio approaches the limit 
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Even for Newtonian liquids, it has not been possible 
to predict the observed diameter ratio from a knowledge 
of the flow velocity and the physical properties of the 
fluid. This is due to the lack of a quantitative expression 

for the axial normal stress at the capillary exit. 

1 

Oliver (68) has measured the thrust exerted by Newtonian 
and non-Newtonian jets and has shown that Newtonian 
liquids also display an axial normal stress under the 
circumstances of a relaxing velocity profile. The usual 
practice has been to calculate this tensile stress from 
the momentum balance and the experimentally determined 
diameter ratio. 

2.2 BREAKUP OP A LAMINAR CAPILLARY JET 

Considerable theoretical and experimental efforts, 
have been made m investigating the laminar breakup of a 
capillary jet of a Newtonian liquid. The major analytical 
contributions have come from the pioneering works of 
Lord Rayleigh (70-73) and Weber (89) . 

Lord Rayleigh (70,71) developed the first under- 
standing of the causes of jet instability through a 

l i 

consideration of pqtential and kinetic energies of an 
mviscid jet moving through a non-interacting medium. He 
showed that if the wavelength of the propagating disturbance 
exceeded the circumference of the jet, the surface tension 
forces were not at a minimum. This condition is unstable 
and results m the breakup of the coherent liquid column 
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into spherical droplets, each separated by the disturbance 
wavelength. Through a numerical solution of the character- 
istic equation, the growth rate of the most rapidly growing 
disturbance was obtained as 

where cy^^is the growth rate for a jet of radius-a^o, 
surface tension- G" and density- £ . The wavelength- £ 

of the disturbance wave corresponding to the maximum growth 
rate is given by 


S = JOZQ * o < 2 * 5) 

If the breakup length is assumed inversely proportional to 
the growth rate, the intact length of the liquid column-L 
is obtained as 

L/za- = c (2.6) 

where Voc is the average jet velocity and C is an exper- 
imentally determined constant. 

Rayleigh's predictions for the breakup length of a 
Newtonian liquid have been confirmed by many investigators 
(75,80,85). His later attempts (72) to include the effects 
of fluid viscosity and interactions with the surroundings 
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resulted in analytical expressions containing Bessel 
function ratios which were not solvable m any convenient 
form. It was possible to conclude, however, that in the 
limit of very high viscosities, the jet would be completely 
stable . 

In 1931, Weber (89) considered the same equations - 
as did Rayleigh, but was able to show that, when the 
disturbance wavelength was much larger than the jet radius, 
the Bessel functions could be approximated by their leading 
term. With this simplification, the growth rate of the 
fastest growing disturbance may be written as 



The presence of viscosity™ is seen to stabilize the jet. 
Corresponding to this maximum growth rate is the wavelength 
of the propagating disturbance which is expressed as 


'2SKtx> 




I * -^222— 


7* 


( 2 . 8 ) 


The length of the intact column prior to breakup is given 
by 


L „ — £ r C F Weif + 3 y s k os, 

ZCU L 

where We oo is the jet Weber- number. We <*#= f an <j 

Re ^ is the jet Reynolds number-Reoo = AH quantities 


(2.9) 
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m equation (2.9) are based on the jet diameter rather than 
the nozzle diameter. The two are usually different due to 
the velocity profile rearrangement. Weber's predictions 
for low speed breakup have been verified by many investi- 
gators (19,33,34,39,50,86,87) . 

The assumptions of the derivation are as follows: 

1. Upon ejection from the capillary, all stresses 
generated by the previous shear flow relax in a distance 
which may be neglected in compansion to the breakup 
length and the jet then moves with a uniform axial 
velocity-Voo (apart from small perturbations) . 

2. The disturbance propagates as a wave which is 
periodic m the axial direction and exponential in time. 

3. There is no interaction with the ambient air. 

4. The disturbance wave is axisymmetric. 

5. The fluid is incompressible. 

In the following paragraphs, detailed comments will be made 
regarding experimental verification of the first four 
assumptions. 

The relaxation of all stresses due to the 
capillary shear flow is accompanied by changes m the jet 
diameter and is considered complete when the jet diameter 
remains constant with increasing distance from the nozzle 
exit. Section 2.1 has previously detailed the ample 
experimental evidence (32,57) proving that, over the normal 
range of Reynolds numbers encountered, the relaxation 
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length xn a lamxnar, Newtonxan jet may be neglected xn 
comparison to the breakup length. 

Donnelly and Glaberson (16) xmposed audxo frequency 
dxsturbances on the motxon of the jet and confxrmed the 
exponentxal growth rate of the propagatxng wave. 

Fenn and Mxddleman (19) establxshed that the 
surround xn^ exert no xnfluence on the jet for Weber 
numbers less than Weoo = 5.3 but, above thxs, a frxc- 
tional interactxon does occur whxch results xn a 
maxxmum xn the breakup length vs. fluxd velocxty curve. 
The locafcxon of thxs maxxmum xs dependent on the ambxent 
pressure . 

Grant and Middleman (33) found that lamxnar jets 
do propagate dxsturbances as an axisymmetrxc wave, xn 
agreement wxth assumptxon 4. Beyond the maxxmum, trans- 
verse waves appear close to the nozzle but the jet stxll 
persxsts xn bexng destroyed by the growth of an 
axxsymmetrxc wave. The nozzle length dxd not exert any 
xnfluence on the breakup length for Jl /2a 0 ratios greater 
than 26. Thexr breakup length data for lamxnar, 

Newtonxan jets are correlated by the express xon. 

L/2a c = 19.5 ( We Q 2 + 3We 0 /Re 0 ) * (2.10) 

l 

where all dimensxonless groups are based on the nozzle 


dxameter . 
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As previously noted, profile relaxation results in a 
change of the jet diameter from that of the capillary bore. 
Kroesser and Middleman used the relaxed diameter instead of 
the capillary diameter to correlate their Newtonian jet 
breakup data and recalculated the measurements of Grant 
and Middleman. The resulting correlation is given by 

, h 0.87 

L/23^ = 18 ( We M + 3We M Ae w ) (2.11) 

The departure from the unit slope prediction of Weber is 
claimed to be statistically valid. 

Weber's analysis assumes the amplitude of the initial 
disturbance to be constant and predicts a uniform droplet 
diameter after breakup. When the disturbance is initiated 
as a result of naturally occurring noise, the initial 
amplitude is a function of time and this periodicity 
produces a droplet size distribution. Duffie and Marshall 
(17) investigated the droplet size distribution to be 
expected from laminar, low viscosity Newtonian jets. Their 
data fit the relation 

D g = 36 ( 2a o ) 0 - 56 ( R eo )-°- 10 ( 2 . 12 ) 

where Dg, the geometric mean diameter, and a Q are m 
microns. For a Sample of N droplets, the geometric mean 
diameter is defined as 
In D g « X(ln Dj/N 


(2a. 12 a) 
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In a second analysis, Weber modified the boundary 
conditions to include the external pressure of the gas. 

His solution predicts the experimentally observed maximum 
in the L vs. V curve but is only in qualitative agreement 
with the data beyond the maximum (19,33). 

Borodin and Dityakm (6) also attempted to evaluate 
the effect of air interaction by writing the equations of 
motion for both the get and the surroundings. Common 

t 

boundary conditions were imposed which resulted m an 18th 
degree equation with complex coefficients that was not 
evaluated because of its complexity. 

For completeness, the paper of Lee and Spencer (41) 
is recommended for its extensive collection of photographs 
showing the effects of injection pressure, air density and 
orifice roughness. A good introductory review of 

investigations into get stability up to 1955 is found m 

1 

Miesse (63) . The empirical techniques used to correlate 
turbulent atomization, data are summarized by Marshall (46) . 

In recent years, considerable theoretical and 
experimental efforts have been directed towards under- 
standing the ungsual behavior of viscoelastic fluids. 

This research has been spawned by the rapid growth of the 
polymer and polymer related industries. 

A fundamental question concerning viscoelastic 
liquids is whether the ability to store elastic energy 
enhances or retards their hydrodynamic stability m a given 
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flow geometry, relative to a .Newtonian fluid. For non- 
Newtonian liquids, the basis for comparison must be in 
terms of its apparent viscosity at the shear rate of the 

t 

experiment. The available evidence is inconclusive as 
data exist to support both conclusions. However, there is 
no reason to believe that viscoelasticity is only 
stabilizing or only destabilizing. It may depend upon the 
magnitude and the algebraic sign of the elastic para- 

j 

meters comprising the particular fluid. 

Wilcox et. ai. (93) found that droplets of 
viscoelastic fluids resisted breakup m a high velocity 
air stream compared to Newtonian liquids of the same zero 
shear viscosity. Davies (12) also noted a stabilizing 
effect by viscoelastic additives on the atomization of 
fluids. 

Dombrowsky, Eisenklam and Fraser (15) studied the 
behavior of a highly elastic Napalm solution m an 
expanding sheet formed by a fan jet. Instead of breaking 
up into droplets as Newtonian liquids do, the Napalm 
sheet formed a large number of filaments. 

Many investigators have studied the onset of 
secondary motions between rotating concentric cylinders as 
their relative velocity is increased for both Newtonian 
and viscoelastic liquids. The instability appears as a 
series of rolls whose onset can be detected visually by 
the use of dyes or mechanically through torque measurements. 
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A review of the extensive results m this geometry is 
considered useful towards appreciating the behavior of 
viscoelastic jets m the capillary jet geometry. 

In 1923, G.I. Taylor (83) presented a linearized 
analysis for a Newtonian fluid contained between concen- 
tric, rotating cylinders and showed that instabilities 
occur whenever the critical value of a certain parameter, 
since named m his honor, was exceeded. His experimental 
results were m excellent agreement with the theoretical 
predictions. Most of the available studies with Newtonian 

I 

liquids for the onset of secondary motions are summarized 
m the books of Lm and chandresakar . When the 
disturbances become finite, the problem must be treated on 
the basis of non-lmear stability theory. 

Many analytical treatments of rotational stability 
with viscoelastic fluids are available using linear and 
non-lmear constitutive relations. Ginn and Denn (28) 
cite or review several of these models m their recent 
paper. The results show that the critical Taylor number 
is dependent upon the magnitude and, m the case of some 

t 

models, the algebraic sign of the fluid's elastic para- 
meters. Thus, the question of whether viscoelasticity 
enhances or retards the hydrodynamic stability, relative 
to a Newtonian fluid of the same apparent viscosity, 
cannot be answered inclusively but is specific to the 
particular viscoelastic fluid. These conclusions are, of 
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course, dependent upon the ability of the constitutive 
equation to describe the real behavior of viscoelastic 
fluids . 

Merrill, Mickley and Ram (49) investigated the 
rotational stability of polystyrene, polyisobutylene and 
polymethylmethacrylate dissolved m various organic 
solvents and found little deviation from Newtonian results 
for liquids of the same apparent viscosity. However, 
for two samples of ultrahigh molecular weight poly- 
lsobutylene, significantly lower Taylor numbers were obtained. 
Rubin and Elata (74) reported substantial stabilization 
for solutions of polyethylene oxide, polyacrylamide and a 
natural polysaccharide. However, viscosity data presented 
for the non-Newtonian solutions is at variance with many 
published measurements. Low zero shear viscosities are 
reported, on the order of 1-6 centipoise for concen- 
trations up to 0.2% and the fluids are claimed to exhibit 
Newtonian behavior except at high shear rates. It is 
suspected that degradation has occurred either m the 
initial preparation or by exposure to the surroundings. 

Denn and Roisman (13) used polyethylene oxide, 
polyacrylamide, hydroxyethylcellulose, carboxymethylcellulose 
and polyisobutylene solutions m their investigation of 
rotational stability. The experimentally determined Taylor 
numbers showed a general trend towards flow stabilization, 
particularly m the narrowest annular gap. In many cases 
though, there was little difference from Newtonian values 
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and there was some evidence of destabilization for poly- 
isobutylene m the widest gap. 

The results of the various experimental investigations 
just described have shown that for selected viscoelastic 
fluids some stabilization or destabilization may occur, 
but, for most of the liquids employed, little deviation 
from Newtonian measurements were found. It is difficult 
to obtain a quantitative comparison between experimental 
data and analytical predictions because the elastic para- 
meters m the viscoelastic constitutive equation^ cannot 
be reliably determined for the dilute solutions of these 
investigations . 

Only limited attention has been given to the breakup 
of a capillary jet of a viscoelastic liquid. Middleman 
(58) presented an analysis for a viscoelastic fluid whose 
behavior is described by a linearized, three constant 
Oldroyd model. The constitutive equation for this material 
is given by 


( 1 + A, ( I + /ft 1^0 A'^ (2.13) 

where Xi and ^M-\ are the relaxation times and retardation 
times, respectively, with X t /!-, . Any variable 

f(r,z,t) was assumed to have the general form 

lkz + °< t 


f (r,z, t) = F(r)e 


(2.14) 
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where k, the wave number, is defined as k=2Tl/S 

i 

Equation (2.13) with the aid of equation (2.14) can be 
written as 


* (2 - 15) 

This result can be regarded as the constitutive equation 
for a Newtonian fluid having the viscosity 


fj -fj- * (2.16) 

The problem is then solved m a manner paralleling the 
development of Weber. 

When X ( and are very small. Middleman showed that an 
expression can be obtained for the breakup length of a 
viscoelastic jet-L v compared to the breakup length of a 
Newtonian jet of the same zero shear viscosity- "fy, . 



Mi 

Q,eo 


Vco • 


1 

+ 3U]eoo/f(eoo 


(2.17) 


Since / a capillary jet of a viscoelastic fluid 

obeying equation (2.13) is less stable than a Newtonian jet 
of the same zero shear viscosity. 

In a recent paper, Kroesser and Middleman have confirmed 
the predicted decreased stability for capillary jets of 
polyisobutylene m tetralm. Quantitative comparison with 
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the prediction of equation (2.17) could not be made without 
known values for A, and Ja, s 

The analytical result of Middleman assumes that the 
jet relaxes from all stresses generated during the previous 
capillary flow m a distance which is short compared to the 
breakup length and therefore regains its zero shear viscos- 
ity 'instantaneously'. He then argues that the shear rates 
associated with the growth of symmetrical disturbances are 
sufficiently low that departure from Newtonian behavior 
would be slight. Thus, the jet propagates disturbances at 
a rate corresponding to its zero shear viscosity. 

As discussed m section 2.1, there have been no 
published measurements on the rate of stress relaxation m 
viscoelastic jets. The data of Gill and Gavis as well as 
that of Middleman and Rautenbach, as reported by Middleman, 
indicate that a slowly decaying axial tension exists 
throughout the lifetime of the jet. Such behavior would 
invalidate an analysis based on an initially relaxed 
viscoelastic column. 

The whole area of stress relaxation and viscosity 
buildup m a newly emergent, non-Newtonian jet will be 
examined m the light of the experimental data taken m 

this investigation. Since the lifetime of a jet is of the 

-2 

order <of 10 seconds, P jet experiment affords the opportu- 

nity to witness effects whose existence is of a similar 



47 

order of .^magnitude. Such short time scales are difficult 
to determine on conventional instruments . 
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CHAPTER 3 

EXPERIMENTAL APPARATUS AND PROCEDURES 

3.1 OVERALL DESCRIPTION OF EXPERIMENTS 

This chapter provides a detailed description of the 
experimental apparatus and techniques. When necessary, 
additional amplification of the procedures will be given 
m later chapters along with the test results. 

A schematic diagram of the apparatus is shown m 
figure 3-1. Either an Instron mechanical tester or 
pressurized nitrogen activated the piston of a hydraulic 
cylinder, which contained the test fluid. The liquid flowed 
through a section of flexible, reinforced hose and past a 
50-125 micron filter. These served to minimize mechanical 
vibrations from being transmitted along the walls of the 
stainless steel piping or through the fluid itself and 
also to trap any contaminant particles. The fluid then 
entered a cube designed to house a pressure transducer and 
the capillary. The cube was mounted on a .concrete block 
to dampen vibrations from the surroundings. 

Still photographs and high speed motion pictures were 
taken of the capillary at magnifications ranging from 
0.8-25 times the actual size. The jet motion was effect- 
ively stopped by the use of microsecond duration flash 
equipment . 

Typical experimental measurements included the jet 
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fnstron 

crosshead 



reservoir 


Figure 3-(. Experimental apparatus for 
studying jet behavior. 
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breakup length, droplet diameter, disturbance wavelength, 
volumetric flow rate and pressure drop through the 
capillary. 

3.2 MECHANICAL EQUIPMENT 

The Instron mechanical tester is a constant speed 

device which moves a crosshead at a preset rate and 

automatically adjusts the force required to do so, up to 

a maximum of 10,000 lbs. of force for the model used. 

Thus, with the Instron crosshead acting on the piston of 

a 4 inch diameter cylinder, containing the test fluid, a 

2 

maximum pressure of 800 lbs/m is available. This type 
of instrument is particularly useful m experiments requir- 
ing a constant volumetric flow rate m capillaries of 
varying length. It was used in capillary viscometer 
experiments and m preliminary tests to determine the 
effect of capillary length to diameter ratio on the jet 
breakup length. 

When pressurized nitrogen served as the prime mover, 
Hoke gas regulators were used to reduce the pressure to 
the desired working value. 

Cylinders of differing construction were used to 
provide low and high pressures. The low pressure set, 
manufactured by the Parker-Hannafm Co., had a 4 inch 
diameter piston by 10 inch stroke. Limitations on 
leakage past the Buna-N piston seals set the maximum 
operating pressure on the steel pneumatic cylinder at 
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250 lbs/in . The hydraulic cylinder was fabricated from 

chrome plated brass and used Viton rubber as the piston 

seal material. It was found that the chrome plating was 

not durable and the exposed base metal corroded if kept 

in contact with the test fluid for several hours. This 

did not pose any ! problem and only required an initial 

flushing with fresh test fluid at the beginning of a 

day's run to cleanse the system. The hydraulic cylinder 

2 

was rated at 500 lbs/in . For high pressure operations 
/ 2 

up to 3000 lbs/m , steel cylinders of 5 inch bore by 12 
inch stroke were used. These were fabricated by the W.E. 
Hennells Co. and employed polyurethane U-cups as piston 
seals m the pneumatic cylinder and Viton rubber f;or the 
same purpose m the hydraulic cylinder. Bronze plating 
was used on all surfaces m contact with the test fluid. 
This lowered the corrosion rate compared to chrome plated 
brass but did not eliminate it. 

Special consideration was given to safely mounting 
the cylinders since the axial force generated on a 5 inch 
diameter piston by a 2000 psig pressure is about 39,000 
lbs. of force. The cylinders were fastened to appropri- 

* c 

ately sized st,eel 'I' beams with high tensile strength 

! 

i L 

bolts. In selecting the components, careful attention 

i i t 

was given to the shearing stress and bending movements 
acting under the most severe loading conditions. Plates 
were welded on either side of the cylinder legs to prevent 
shifting during any cycling operation. 
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A Dynisco 0-1000 psig pressure transducer was used 
to measure the pressure drop of the fluid flowing through 
the capillary. The millivolt output of the transducer was 
plotted on a Leeds and Northrup 'Azar' recorder. 

Nozzles were fashioned from hypodermic needles or 
tubing by grinding their ends to a flat face. The bores 

i 

were examined under 100 power magnification for 
concentricity and surface imperfections. A stage micro- 
meter was used to determine the inside diairieter. The 
length/diameter ratio of the nozzles were usually greater 
than 100/1 although, m certain experimental runs, shorter 
nozzles were used to 'maintain a safe working pressure. 

3.3 PHOTOGRAPHIC EQUIPMENT 

Still photographs were taken with either of two 
4"x5", triple extension (22") bellows cameras made by the 
Bausch & Lomb Co. and the Calumet Photographic Co. A 

} r 

32 mm. f4.‘5 lens was used for closeup work of high 
magnification and a 158 mm. f6.3 lens was used for general 
purpose photography of lower magnification. The cameras 
were mounted on a specially constructed table which allowed 
movements along the x,y and z axes. 

Both 3''x4" Polaroid roll film and ^"xS" Polaroid 
sheet film of ASA 3000 rating were used. The high speed 
film allowed small lens openings to be used with the 
accompanying advantage of a wide depth of field. 

Motion pibtures were made with a 400 foot capacity. 
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16 miru'Hycam 1 series camera manufactured by Red Lake Labora- 
tories. It is capable of frame rates between 10-11,000 

t 

pictures per second. By means of a timing light, which 
placed a mark on the edge of the movie film at a 
frequency of 10,100 or 1000 pulses per second, the 
instantaneous frame rate could be determined. An 
accessory event sychromzer fired an electronic flash 
when the 1/2.5 shutter was fully open. 

The microsecond duration, high intensity lighting was 

i 

provided by electronic single flash and strobe units pror~ 
duced by Edgerton, Germeshausen and Grier, Inc., models 
549 and 501 respectively. A General Radio Co. strobe, 
model 1538, was available for visual measurements under 
circumstances wljere photographic determination was 
difficult. 

3.4 GENERAL EXPERIMENTAL PROCEDURES 

The experimental procedures used to obtain the test 
results are generally straightforward and do not require 
extensive explanations. An overall description of the 
techniques that were utilized has .been given previously in 
section 3.1- . Some additional details will be furnished in 
this and the following sections. 

When nitrogen gas was used' as the prime mover instead 
of the Instron, two to four holding tanks were first 
pressurized to the desired working level. The gas regulator 
was then valved off from the system in order to avoid the 
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possibility of periodic disturbances produced by the 

s' 

opening and closing of the regulator. The decrease in 
pressure during an experimental run, due to piston travel, 
was negligible. 

A Dynisco strain gage transducer measured the fluid 
pressure at the entrance t:o the capillary. The 
calibration wars checked periodically against a dead weight 
tester and found to be invariant. A typical plot of the 
transducer's millivolt output versus applied pressure is 
shown in figure 3-2. 

Volumetric flow rates of the test fluid were 
determined easily with a graduate cylinder and a stop- 
watch. When the Instron was used in the experimental 
system, the crosshead travel could be read to the nearest 
0.005 inches. This movement was translated into volumetric 
flow rates from a knowledge of the hydraulic cylinder and 
capillary diameters. 

3 16 PROCEDURES FOR PROFILE DECAY EXPERIMENTS ! 

t 

Accurate measurements were taken of the 'jet 
diameter as a function of the axial distance in the region 
immediately outside of the capillary. This is related to 
the rate at which the velocity profile decays. 

i t 

The electronic flash was placed behind the jet with 
its light output directed onto a translucent plastic 
disk. A shadow photograph of the capillary ^et was obtained 
at a magnification of approximately 10/1. In order to 
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to correct for any distortion from the camera lens or arising 
from a lack of parallelism between the plane of the film and 
that of the pet, a rod of known diameter was also photo- 
graphed at each location* The required measurements were 
performed on a Nikon profile projector, model 6C, fitted 

l 

with 'a 10/1 enlarging lens and a micrometer stage readable 
to 0.0001 inches. Values reported for the jet diameter 
are probably accurate to more than three significant 
figures. 

Below a minimum velocity, the pet was either 
' ,* 
affected by gravity or the fluid wetted the capillary 

I 

«• 

face to form a larger diameter than the bore. These runs 
were rejected. The temperature of the jet was measured 
after collection and compared to ambient for evidence of 

viscous heating.' A Brookfield viscometer was used to check 

\ 

* 

the viscosity of Newtonian test fluids after collection. 

The shear dependent viscosity of the non-Newtonian liquids. 

i * 

made a similar check on them impractical. 

3.6 PROCEDURES FOR JET BREAKUP EXPERIMENTS 

Still photographs were taken of the capillary jet 
every few centimeters along the liquid column, starting 
from the nozzle tip and ending after the appearance of 

i 

free droplets. A white reflector card was placed 6^12 

* 

inches behind the jet and the area of interest illuminated 
by the EG&G electronic flash, placed at a 45° angle to the 
camera. Figure 3-3 illustrates the^photographic arrangement. 
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Still pictures were taken with the room darkened and the 
shutter speed set on 'time'. A lens opening of f32 was 
used for magnifications of 3/1 while at enlargements of 
25/1, a lens opening of £4.5 was required. 


White Reflector Card 


Nozzle 


Capillary Jet 



Camera 



Electronic Flash 


Figure 3-3. Photography of a capillary jet. 

By means of a plumb bob suspended from the camera to 
a millimeter ruler, the distance between individual photo- 
graphs could be determined, as illustrated m figure 3-4. 

Photograph 1 Photograph 2 




Nozzle k 


Jet 



Ruler 


Figure 3-4. Measurement of distances between photographs. 
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The breakup length of a capillary jet is the 
distance from the nozzle tip to the point where the liquid 
column is no longer coherent and forms individual droplets. 

To determine this, quantity from photographs taken at the 
two locations first requires obtaining the overall distance, 
as described m the previous paragraph. This length is 
then adjusted for the increment that the nozzle tip and the 
breakup point are each offset from the midpoints of their 
respective photographs. A photograph of a millimeter 
ruler placed m the same plane as the jet allowed actual 
distances to be measured on the experimental photographs. 
Ordinarily, an enlargement of 3-6 times actual size was used. 

Moderately concentrated solutions of the visco- 
elastic additives formed capillary jets which often had 
breakup lengths larger than the 5-6 foot span that could be 
conveniently photographed. In these instances, visual 
measurements were made with the aid of the General Radio 
strobe . 

In order to determine the disturbance wavelengths 
and the droplet diameters, more accurate techniques were 
required. The photographs of these subjects were placed m 
an overhead projector and the resulting image cast upon a 
screen had bn overall magnification of 10-11 times actual 
size. Lengths were read off the screen with the aid of a 
dial caliper to an accuracy of 0.001 inches. 
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CHAPTER 4 
TEST FLUIDS 

4.1 SELECTION OF TEST FLUIDS 

The test fluids used in this research can be broadly 
classified into three rheological categories: 

1. Newtonian 

2. Non-Newtonian inelastic 

3. Viscoelastic 

Newtonian fluids consisted of water, ethylene 
glycol and glycerin-water solutions. These were included 
m order to verify the soundness of the experimental 
apparatus and procedures by comparing the data with pub- 
lished results of other investigators. 

The non-Newtonian inelastic fluids were comprised 
of Carbopol 934-water solutions and dispersions of the 
following materials- MPA 60 m xylene, MPA 60 in odorless 
mineral spirits and silica EH-5 in water. Despite their 
compositional differences, the liquids shared certain 
similarities from a rheological point of view. All 
exhibited a yield stress, a shear dependent viscosity and 
an absence of viscoelastic phenomena, e .g. a Weissenberg 
effect or recoil upon removal of shearing stresses. 

Carbopol is a carboxypolymethylene polymer, thought 
to be slightly crosslinked (21) , and is one of the few 
organic polymers whose solutions do not display elastic 
effects. Huppler (37) and Hurd (38) have characterized 
Carbopol solutions as a non-Newtonian fluid which possesses 
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very small normal forces. Several investigators, including 
Dodge and Metzner (14) and Kapoor (40) have demonstrated 
that Carbopol solutions are not drag reducing nor do they 
demonstrate any recoil upon removal of shear stress. 

MPA 60 is a hydrogenated castor oil derivative which 
was selected m order to observe the behavior of a 
capillary jet where the fluid's gellation time clearly 
exceeded- the jet breakup time. They are time dependent 
fluids whose rheological behavior depends upon the previous 
shear history. Upon cessation of shear, the consistency of 
the fluid increases with time and attains a steady state 
structure after about a day's standing. On the other hand, 
the Carbopol solutions and the silica dispersion appeared 
» to gell immediately after the shearing stress was discont- 
inued . 

The aqueous viscoelastic solutions were prepared from 
a class of polymeric additives which exhibit drag reducing 
properties (23,36,76,84) m turbulent flow. These test 
fluids consisted of Separan AP30- a partially hydrolyzed 
polyacrylamide. Polyox: coagulant grade- a polyethylene 
oxide polymer, SCMC 7HS- a high viscosity grade of sodium 
carboxymethylcellulose and Guar Gum- a naturally occurring 
polysaccharide . 

4.2 PREPARATION OF TEST FLUIDS 

All aqueous polymer solutions were prepared by dissol- 
ving the additive m distilled water using one minute of 
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mixing with a laboratory mixer followed by gentle rolling 
overnight. This procedure was chosen to minimize shear and 
thermal degradation. Exposure to outside ultraviolet 
radiation was also avoided as many organic polymers are 
adversely affected by such treatment. The MPA 60 dispersions 
were used as supplied by the manufacturer. Silica dis- 
persions require high shear mixing to break up the 
agglomerates and this was accomplished by batchwise 
addition to a Waring blendor. A non— ionic surface active 
agent, Triton X-100, was added to improve the wettability 
of the silica particles. 

Several of the organic polymers required the use of 
additional additives m order to minimize chemical degrad- 
ation. In private communications. Dr. B.D. Marsh advised 
the use of small quantities of glycerin with Separan 
solutions to offset the effects of carbon dioxide absorption. 
The manufacturer of Polyox recommends isopropyl alcohol m 
solutions to prevent chain sission reactions. Formaldehyde 
was added to SCMC and Guar Gum solutions to avoid biological 
attack. 

4.3 PHYSICAL PROPERTIES OF TEST FLUIDS 

Surface tensions of non-Newtonian fluids were deter- 
mined on a Rossano tensiometer while those of Newtonian 
fluids were taken from published data. The common static 
techniques were found to be inadequate for fluids with 
appreciable yield stresses or those of high viscosity as 
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the measuring element could not penetrate the liquid sur- 
face. For such solutions, a dynamic method taken under flow 
conditions would appear to be necessary. Surface tension 
values for the two MPA 60 dispersions were determined by 
heating a sample to a temperature of approximately 70— 90°C, 
at which point they were substantially mviscid, and 
extrapolating the data back to room temperature by 
assuming the same linear temperature dependence as for 
benzene and toluene. The results are tabulated in table 
4-1. This technique was ineffective with 0.6% Carborpol or 
With the silica dispersion, both of which remained gellified 
at elevated temperatures , 

Table 4-1. Surface tension of MPA 60 dispersions 


as a 

function of temperature. 

Fluid 

Temperature 

Surface tension 

MPA 60 in xylene 

80 C 

24.6 dynes/cm. 


75 

25.6 


70 

27.2 


25 (est.) 

32.5 

MPA 60 in 

90 

19.0 

odorless mineral 

85 

19.4 

spirits 

80 

19.7 


75 , 

20.2 


25 (est.) 

26.0 


Surface tension measurements on the moderately 
concentrated viscoelastic solutions were found to depend upon 
the rate at which the platinum blade was withdrawn from the 

5 

solution, with slower rates producing lower values. This 
variation is probably related to relaxation phenomena. 


Within the same concentration range as that used here, 
the density of aqueous solutions of the organic polymers has 
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been shown by many authors to be equivalent to pure water 
and this value is accepted . 

Table 4-2 is a classification of the test fluids 
and lists some pertinent physical properties. All 
Newtonian properties were obtained from published values 
m the standard handbooks. For the convenience of the 

i 

reader, the density, surface tension and viscosity of 
glycerin-water solutions are shown graphically as a 
function of composition m figure 4-1. 

4.4 DEVELOPMENT OF THE COAXIAL CYLINDER EQUATIONS 

A coaxial cylinder viscometer was used to obtain 
shear stress-shear rate data at Low to medium rates of shear 
and was also employed m the determination of yield 
stresses. The equations of flow for this geometry will 
be derived in a manner similar to that presented by 
Middleman (59) . A schematic representation of this 
instrument is shown m figure 4-2. 

A torque-M rotates themner cylinder of radius-R^ and 
height-H at an angular velocity-./! . The shear stress at 
the wall of the bob is given by 

_ M 

b 2 TTr£(H+H 0 ) (4,1) 

where H 0 is the end correction of the rotor. The outer 
stationary cup of radius-R c is under a shearing stress- 
which is written as 



Table 4-2. Physical 


Test fluid 

Carbopol 934 
( carboxypoly- 
methylene :B ,F . 
Goodrich Co.) 

Ethylene glycol 
Glycerin-water 

Guar Gum 
(natural poly- 
saccharide: 
Stem-Hall Co.) 

MPA 60 

(hydrogenated 
castor oil deri- 
vative: Baker 
Castor Oil Co.) 

SCMC 7HS 
(sodium carboxy- 
methylcellulose: 
Hercules Powder 
Co) 

Separan AP30 
"(hydrolyzed poly- 
acrylamide: Dovt. 
Chemical Co.) 


Concentration 

cl 

0.1% m water 

' 3 , 

0 1 2% m water 

a 

0.6% m water 
100 % 
varied 

0.05% in water‘d 

24% in xylene 

24% m odorless 
mineral spirits 

c 

0.25% m water 

G 

0.7% m water 

0.05% m water^ 
0 . 25% in water 


Densit\ 


1. gm./cc. 


1.116 

see figure 
4-1 


0.872 


0.816 


1 . 


of the test fluids 



32.5 dynes/cm. 
26.0 

72 

72 


see table 4-3 
see table 4-3 

0.45 poise 
2.3 


72 0.50* 

CFi 

63-72 (time 10.0* 

dependent) 



Table 4-2. Physical properties of test fluids (continued) . 


Test fluid 

Concentration 

Density 

Surface tension 

Newtonian or zero 
shear viscosity 

Silica EH-5 
(Cabot Co.) 

0 

12.4% m water 

1.16 gm/cc 

not measurable 

see table 

4-3 

Polyox-coagulant 
(polyethylene 
oxide: Union 
Carbide Co.) 

f 

0 . 25% m water 

1 . 


2-5 poise 

(est. ) 

Water 

100% 

0.997 

72 dynes/cm. 

0.0092 



a- neutralized with 0.4 gm. of sodium hydroxide/gm. of Carbopol. 
b,c- formaldehyde added to prevent bacterial growth, 
d- small quantities of glycerin added for stability, 
e- Triton X-100 (Rohm and Haas Co.) wetting agent added, 
f- isopropyl alcohol added for stability, 
data of Schwarz and Bruce (77). 


U1 
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Figure 4-2. Typical coaxial cylinder viscometer, 


T c 




ETTR^Ch + M*) 


(4.2) 


If the fluid has a yield stress- , there are three 
possible, conditions of flow behavior. 


“Cu 4 T c 


(4.3a) 

(4.3b) 

(4.3c) 


In the condition demanded by (4.3a), there will be no flow 


m the gap. When inequality (4.3b) holds, the fluid will 
flow only m that portion of the annular space where the . 
local shear stress exceeds the yield stress and the effective 
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gap size will be less than that defined by the viscometer 
walls. The critical radius at which flow occurs is obtained 
from the relation 

M-zITR^CH + HoV^ ( 4 * 4 ) 

which is rearranged to 

^ c (4 ’ 5) 

By virtue of equation (4„1), the critical radius for flow 
may be written as 


- R, fE 

The minimum shearing stress needed to insure flow across 
the entire annular space will occur when . From 

equations (4.1) and (4.2) 


Tc 


,1A\A 


or (Cb) 


VAvA. 


ClTb\y uv>\ — Ssi 

Xj Rb 



(4.7) 

(4.8) 


In a simple shear flow field where the only non-zero 
velocity component is V© = VeC r ) ■= X \aJ ( r) , the rate of deformat 
ion tensor, Av^ , is given by 


^ er 




(4.9) 


where uJ (r) is the angular velocity. 
From the equations of motion 
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r 't'er = constant 


(4.10) 


Taking the total differential of equation (4.10) and 
dropping the subscripts 


2 rtdr + d 't «o 


(4.11) 


or 1= - 21 A- 

dt r dtr 


(4.11a) 


Applying this transformation to equation (4.9) 

£cp = -i - r ^ = 


(4.12) 


Integrating and applying the boundary conditions that uj =o 


at r=R c and that W =SL at r=Rj- ) 


St - 






(4.13) 


Equation (4.13) is differentiated with respect to to yield 


-da = _l itci 

dX b 2 [ th ^ dr b j 


(4.14) 


This result is valid only if the shearing stress at the cup 
exceeds the yield stress of the fluid 


From equations (4.1) and (4.2) 
Therefore 

dr b \Rc ; 

4 

and equation (4.14) becomes 


(4.15) 


(4.16) 


zx b = -?cc b ) - p 


(4.17) 



70 


If R C ^R^ or, equxvalently, xf s<3(l a closed form solution 
of equatxon (4.7) xs obtaxned as 


- ***&. = 25t < 4 - 18 > 

The general case for fxnxte s is treated by an xterative 
procedure. Equatxon. (4.18) is valxd for all values of , 
e.g. etc-, and may therefore he wrxtten as 


4ik 

dr b 

z 'tb 4 ^- 




<ix b J s n^ b 


z. ^ 4^- 


= -CtV-tto) - -r(S H Xb) 


(4.19a) 

(4.19b) 


(4.19c) 


Addxng equatxons (4.18) and (4.19a) -(4.19c) produces the 
summation equation 

* <■*■> - m t ** (4 - 20) 

As N -^»oo , -f-CS 2- ^) — ^ Q 


. The result obtained is 


C*> 


•Kro = 2t s r s" ^ 


i-o 




or 


^■CTv^ — 2LJL V~ s iL AlsaJk 


(4.21) 

(4.22) 




If g- d In it xs constant, as it xs for a power law fluid, 
d In Xio 

xt may be shown that the shear rate at the bob xs gxven by 


t W; 


2.SL T~ 

C— O 

f — 




% 






(4.23) 

(4.23a) 
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4.5 COAXIAL CYLINDER MEASUREMENTS 

Shear stress-shear rate c3ata at low to medium rates 
of shear and yield stress measurements were obtained on 
the Stormer and Brookfield LVT viscometers, modified so as 
to provide mathematically treatable results. The Brookfield 
guard ring was removed and simple adaptions on both instru- 
ments permitted the use of the rotors and cup assembly 
from the MacMichael| viscometer. 

End corrections for the rotors were established with 
a Dow Corning 50 cs. silicone oil, which has been shown to 
be Newtonian up to shear rates of 10^ sec”'*' (19) . The 
method outlined by Van Wazer et.al, (88) was used and involves 
measuring the angular velocity- as a function of the 
torque applied to the rotor at several values of H, the 
depth of immersion. Experimental data are plotted as M/f|_ 
vs. H, which is a straight line for Newtonian liquids. The 
line is extrapolated back to (m/J"L ) =0 and the negative 
intercept-H 0 corresponding to this point gives the corrected 
height of the rotor as (H4H 0 ). Figure 4—3 is a typical plot 
of the data taken on a modified Stormer viscometer for a 
2 cm. diameter rotor in a 3.5 cm. cup and the extrapolat- 
ion shows H 0 to be 0.38 cm. 

The yield stress for a test fluid was determined on 
the Stormer viscometer by noting the minimum weight , required 

to turn the rotor. From the dimensions of the instrument, 

' 2 
the yield stress in dynes/cm 


is given as 
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Figure 4-3. Determination of the end correction for 
a rotor in a. modified 5tormer 
viscometer; 

test fluid- 50c.s. silicone oil, bob 
radius- I cm 5 cup radius- 3.5 cm. 
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1.425X980 

n- . — - 11 w 20.2 W m;Ln 

C'j - * w mm = 5 (4.24) 

2 r£ (H+H q ) R b (H+H 0 ) 

where is the minimum weight which will turn the rotor. 

Yield stresses were also obtained on the Brookfield visco- 
meter by twisting the rotor until the dial reading is at 
full scale and then releasing it. The dial reading 
where the pointer ultimately stops may be converted to 
a yield stress by the relation 

^ fraction of full dyne-cm torque for 

'T!\j = ecale reading full scale reading (4.25) 

2 TrR b (H+H 0 ) 

A final steady state value for the yield stress was consi- 
dered to be reached when the instrument reading remained 

f 

constant for several hours. This criteria will not detect 
any fractional revolution of the rotor at rates slower 
than this time period. 

Table 4-3 lists the yield stress determination 
for the viscoinelastic test fluids. Fisher et.al. (21) have 
also measured yield values for Carbopol solutions and their 
results are plotted along with the current value in 
figure 4-4. Good agreement is seen between the data points 
of Fisher and that of this investigation. Since the MPA 60 
dispersions have long gelation times, they were prepared 
for testing by pouring a freshly sheared sample into the 



YIELD STRESS - dynes /cm 
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8 


CARB0P0L CONCENTRATION- weight % 

Fi'0u.re4-4. Yield stresses of Carbopol solutions; 

O-dafa of Fisher etal. (“21 ) ? A- the 
present worK. 
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annular space between the aup and rotor, then allowing 

9 

the thixotropic fluid to set for 2-3 days* An appropriate 
cover was provided to prevent solvent evaporation. The 
resultant accuracy of the yield stress for these fluids 
is uncertain because of the tendency of the MPA 60 to 
separate from the rotor while undergoing gellation. 


Table 4-3. Yield stresses of the viscomelastic 
test liquids. 


_ Test liquid 

0*1% Carbopol 
0.2% Carbopol 
0.6% Carbopol 
24% MPA 60 m xylene 
24% MPA 60 in mineral 
spirits 

12.4% silica m water 


Yield stress 
0 dynes/cm^ 

55 

800 

800 


2000 


Shear stress vs. angular velocity measurements at the 

I 

rotor were obtained on the Brookfield and S former viscometers 
for 0.1%, 0.2% and 0.6% Carbopol solutions. Similar data 
could not be acquired with the MPA 60 dispersions because 
of their thixotropic nature and is also absent for the 
12.4% silica dispersion as this highly gellified fluid 
separated from the rotor, once motion was initiated, and 
did not flow back. The measurements were plotted m figure 
4-5 and show that the fluids exhibit power law behavior. In 
all cases, the data were examined to insure that the shearing 
stress at the cup exceeded the yield value, a procedure 
which guarantees flow across the entire gap. The slope of 
the lines, n= d log 7^/ d log -ft, , is given m table 4-4. 


SHEAR STRESS AT ROTOR- Z b dynes/em 



ANGULAR VELOCITY OF ROTOR - n md,/sec. 


Figure 4-5, Coaxial cylinder measurements in modified. Sformer and. BrooKfielct 
viscometers; x~O.G>%Carbo<po\ ) 0-0.1% Carbopol, n-0.\°/o Carbopol *, 
continuedon next page. 




SHEAR STRESS AT ROTOR- X b dynes/cm 


JO 1 



ANGULAR VELOCITY OF ROTOR - S 2 . md./sec. 


Figure 4-- 5. corf mued. 


-j 

-4 
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Table 4-4. Power law exponent for Carbopol 
solutions from coaxial cylinder 
viscometer measurements. 


Fluid n 
0.1% Carbopol 0.88 
0.2% Carbopol 0.32 
0.6% Carbopol 0.20 


Figure 4-5 results will be combined with the high 
shear rate data from the capillary viscometer m section 
4.8 to provide shear stress-shear rate data over a wide range. 


4.6 DEVELOPMENT OF THE CAPILLARY FLOW EQUATIONS 

Viscometric data at high shear rates were obtained 
with a capillary viscometer. 

The classical Rabinowitsch derivation for the wall 
shear rate of a general fluid will be modified to show that 
the existence of a yield stress does not alter the final 
result. A typical section of a capillary tube is shown m 
figure 4-6 „ 


Pressure exerted 

on face = P + A P Pressure 



Figure 4-6. Force balance for a capillary 
viscometer. 
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The pressure- AP acting on the ends of the fluid 
cylinder exerts a force-F equal to 


f = AP *1 Tr 2 - 


(4.26) 


If the fluid possesses a yield stress, viscous flow will 
occur only over that portion of the capillary crossection 
where the shear stress exceeds the yield value; other- 
wise, plug flow results. Above the critical radius for 
flow-ry,-the force arising from the pressure is resisted 
by the viscous restraining forces 

F = tr- ztrrk (4 ' 27) 


where 't is the shearing stress. In steady flow, the 
two forces are equal 


V - <->r>j (4.28) 

At the capillary wall, equation (4.28) may be written as 


r w £*(4.29) 

Dividing equation (4.28) by equation (4.29) gives the 
linear relation between the shearing stress and the radius. 


X _ v" 

r* ~ <x 0 


r > ^ 


(4.30) 
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As indicated, it is valid for all values of the radius 
greater than the critical radius. 

The flow rate through the capillary tube is given 
by 


<3 = f o a ° 2Trr VzrCOti'f' 


(4.31) 


Integrating by parts and invoking the zero slip condition 
atr r= a Q 


9 = a^) d,r 

or q - C irr 1 (- ^)dr + 5 ^ IT r 2- (-|S£) dr 


(4.32) 

(4.33) 


But = ° f ° r ° < r < ^ 


\ 

and equation 


(4.33) reduces to 


q = S^irrHCt)^ (4.34) 

where f(X) has been defined as 

<A<f~ 

By virtue of equation (4.30), (4.34) may be written as 


4 Q 

■wot 




^(X)^dX 


(4.35) 

V 


Differentiating equation (4.35) with respect to X 


w 



The last term on the right hand side of (4.36) is identi- 

✓ * 

cally zero since f (^y )=0. Rearranging equation (4.36) 



with, the aid of (4.35) results in an expression for the wall 
shear rate for a general fluid m terms of macros copically 

i 

measurable variables. 


A. = _ 

1 


(4.37) 

Viral/ 


OTw 

or 



(4.38) 

where q is defined as q 

= d ln(4Q/TT a 0 3 ) 

(4-.39) 



d InOT 

w 


It is usual to 

plotT w vs 

. ^ Q m . and to 

iral 

define the slope 

as n = 1/q 

. Equation (4.38) is 

alternately 

written as 




* CM = 



(4.40) 


Many pseudoplastic fluids are found to obey the relation: 


'-u w - 


Apgp 


- K ' /-JU_ \ n 


(4.41) 


This may be put into the standard form for the power law, 


Vw = K - \<X 

by the relation k= K' ( 

\ 3 h + 1 


(4.42) 


(4.43) 


The entrance and exit losses for flow through a 

« 


capillary may be eliminated, by taking pressure drop 
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measurements wxth capxllarxes of two dxfferent lengths but 
at the same flow rate. Thxs method, orxgxnally used by 
Couette, is dxscussed xn Fredrickson x (22) . Consxder a 
capxllary of length- -5 L xn whxch the pressure drop xs AP-q 

and a short capxllary xn whxch the correspondxng quantx- 

6 

txes are _? s and AP S , respectxvely . The total pressure 
drop through the capxllary may be wrxtten as the sum of 
an entrance loss-AP ent , an exxt loss - AP ex and a pressure 
drop due to fully developed flow -APf. 


A P T = A P .+ 
L ent 

(AP f ) L 

+ AP 

= Ap ent + 

< Ap f>s 

+ AP 


ex 


(4.44) 

(4.45) 

Since the flow rates are xdentxcal, the entrance and 
exit losses are the same for each capxllary. Subtractxng 
equation (4.46) from (4.44) gxves the corrected pressure 
drop over the length- ) as 


= APl ~ 

A l “ A s 


(4.46) 


4.7 CAPILLARY VISCOMETER MEASUREMENTS 

t 

Shear stress-shear rate measurements were obtaxned 
wxth a capxllary vxscometer. The Couette two capxllary 
method just described was used to elxmxnate entrance and 
exxt effects. i^.n Instron mechanxcal tester was used to 
drxve the fluxd from the hydraulxc storage cylxnder 
through the capxllary. The constant speed feature of thxs 
xnstrument makes xt highly convenxent for thxs type of 
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experiment, where xt xs necessary to achxeve a fxxed 
volumetrxc flow rate in capillarxes of varyxng dxameter 
and length. Pressure drops were measured wxth a calxbrated 
straxn gage transducer whxle flow rates were accurately 
known from the Instron crosshead travel. A more detaxled 
dxscussxon of procedures xs avaxlable xn chapter 3. The 
nozzles were suffxcxently long to ensure fully developed 
flow. 

The orxgxnal test data for the fluxds, taken xn the 

capxllary vxscometer, are tabulated xn append xx A . 

From these measurements, plots were prepared of the wall 

shear stress as a functxon of the apparent shear rate and, 

as seen xn fxgure 4-7, the fluxds exhxbxt power law 

behavxor. Table 4-5 summarxzes the power law constants- K 1 

n 

and n from the relatxon: a Q AP/2 JL = K* (4V 0 /a 0 ) 

Table 4-5 Power law constants from capxllary 


vxscometer data 

Test fluxd KJ n 

24% MPA 60 xn 8.9 dynes-sec^ 0.74 

xylene cm 2 

24% MPA 60 xn 3.8 0.79 

mxneral spxrxts 

12.4% sxlxca-water 22.5 0.60 

0.2% Carbopol 11.7 0.59 

0.6% Carbopol 504 0.35 


In sectxon 4.8, the capxllary vxscometer measurements 
and those from the coaxxal cylxnder wxll be recalculated to 
obtaxn shear stress-shear rate curves. 



(a 0 aP/ 2Z)- lbf/ft 



Figure4-7 Capillary viscomeier data for feet fluids; 

0-0.0'2G»3 cm. nozzle, X-0.0414 cm. nozzle, A-0 0808 cm. nozzle, 
□ ~0.'Z / 28 cm. nozzle; 
continued, on next page. 


CO 

4 ^ 



(a 0 ^p/zX)-ib f /n 



Figure 4-7. continued 


CO 

in 
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4.8 SHEAR STRESS-SHEAR RATE DATA 

The true shear rate for data taken with the coaxial 
cylinder viscometer was obtained from equation (4.23a) . 

This requires the arxgular velocity of the rotor- , the 
ratio of the bob radius to the cup radius— s and the 
slope of the double logarithmic plot: log_TLvs. logTr^. The 
shear stress at the rotor is given by equation (4.1) . For 
fluids with a yield stress, the data should be trimmed to 
exclude any measurements yhere 'Csj>'tTc . 

In a similar manner, equation (4.40) gives the true 
shear rate for data taken with the capillary viscometer. 

The input required for this equation is the flow rate-Q, 
the capillary radius- a Q and the slope of the double loga- 
xithmic plot: log / t w vs. log 4V 0 /a 0 .. Equation (4 . 29) is used 
to obtain the- wall shear stress. 

The shear stress-shear ■ rate curves obtained for the 
test fluids are given m figure 4-8. Data for 0.2 % 

Carbopol and 0.6% Carbopol show a 'flattening' of the 
curve at low shear rates as the yield stress becomes a 
significant portion of the total applied stress. At higher 
shear rates, the usual power law behavior is indicated. 

Shear stress-shear rate measurements are available 
m the published literature for most of the viscoelastic 
fluids used in this investigation- 0.05% and 0.25% Separan 
(77,78), 0.25% SCMC '(48), 0.7% SCMC (52) and 0.05% Guar Gum (91). 



SHEAR, ST RESS-Z dynes/cm 
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SHEAR RATE- t sec“ s 

Frqure 4-8. Shear stress- shear rate relations for 
test fltxid-Sj 

^-24% MPA 60 m xylene, ©-2.4% MPA 60 in 
mineral spirits^ -12,4% silica in water, 
0-0.S%Carbopol ? A '-A- 0.2% Carbopol, D 
0,6 %Carbopol, unf \ i led symbols for coaxial 
cylinder data, fined symbol s -for capillary 
viscometer data; 
continued on next page. 



Figure 4-8. confinued. 
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CHAPTER 5 

PROFILE RELAXATION OF A LAMINAR CAPILLARY 
JET OF A NON-NEWTONIAN. LIQUID 

When a laminar capillary jet of liquid emerges from 
the nozzle into ambient air, it experiences a sudden removal 
of the confining wall. As a consequence of this, all 
stresses resulting from the previous shear flow will decay 
with increasing axial distance and the velocity profile 
will tend towards a uniform flow. Whether or not, the 
capillary jet is ever completely relaxed within its life- 
time depends upon the rate of stress relaxation. This m 
turn is a function of the rheological character of the 
fluid and the flow velocity. 

Stress relaxation results m a diameter change which 
may be larger, smaller or unchanged from the nozzle bore. 
Experimentally, a jet diameter which remains constant with 
increasing axial distance is taken as evidence that the 
stress relaxation is complete. 

The mechanism of profile decay in non-Newtonian jets 
is complicated by the fact that these fluids possess a 
three dimensional structure whose spatial configuration 
depends upon the local rate of shear. For example, polymer 
molecules m solution are believed to exist as a randomly 
entangled network of chains bound together by secondary 
chemical bonds. Under shear, these macromolecules tend to 
elongate and orient themselves with their major axes 
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parallel to the flow direction. 

At the high shear rates of capillary flow, the three 
dimensional, 'zero shear' structure of the fluid is largely 
destroyed and a finite time is required outside of the 
capillary before it is regained. If the liquid possesses 
a structural formation time which is short compared to the 
time required for stress relaxation, then profile decay 
should occur at a rate corresponding to the liquid's zero 
shear viscosity. At the other extreme, a fluid with an 
' infinite ' structural formation time will relax at a rate 
corresponding to its apparent viscosity at the exit plane 
of the capillary. 

This chapter investigates profile relaxation in 
capillary jets of non-Newtonian inelastic fluids and in 
viscoelastic liquids. 

5.1 THE RATE OP PROFILE RELAXATION 

The prediction of the rate of profile decay in 
Newtonian ^ets is facilitated by separating the analysis 
into two parts. A boundary layer approach is taken at large 
Reynolds numbers for short axial distances where the 
periphery of the liquid column does not interact with the; 
core (30) . At large axial distances, a linear perturbation 
approach is applied about the final relaxed state (32) . 

Goren and Wronski present experimental results that are in 
fair agreement with their theoretical predictions. 

The present effort extends the theory at large axial 
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distances to a capillary jet of a linear viscoelastic 
fluid. 

5.1-1 SOLUTION OF THE DYNAMIC EQUATIONS 

At some distance ’infinitely far' from the nozzle 
tip, the jet will have completely relaxed and attained 
its final values of diameter, velocity and pressure, as 
shown in figure 5-1. The rate of profile decay at large 
axial distances will be determined by performing a small 

f 

perturbation analysis about the final state in; a manner 
paralleling the development of Goren and wronski. 
Frictional drag from the surroundings and gravity effects 
are assumed negligible. The flow is steady and the decay 
progresses in an axisymmetric pattern. 

✓ 



Figure 5-1 Profile relaxation of a laminar 
capillary jet. 

\ 

In the usual notation, at z= o° 

V r =0 ; V z =Veo ; a=a 

h, ! r if° 


(5.1) 
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Just before z- oo 

V r = V^; V z '= Voo + V^, ; a = a*, + a 
P = Poo + p" ; T^j = 0 + t 


(5.2) 


The perturbation quantities /designated as ( ) / are 
regarded as being infinitesimally small . Thus , quantities 
which are second order m magnitude , such as the inertial 
terms in the equations of motion or^other non-linearities 
m the rheological constitutive equation, can be neglected. 

The linearized equations of continuity and motion 
m a cylindrical coordinate system can be written as 


-Sr £= CrlM 


r ar 

\Ioo 

oV oo ^ — _ lt£L 

v 32 32 


dVr _ 

3 2' 


3Ve 

32 

dr 


O 


_L 3_-(rTrr)_ & 4 . 

r ar ' r 32 
r ^ a? 


(5.3) 
(5 . 4a) 
(5.4b) 


The boundary conditions state that the shear stress at \ 
the }et surface is zero since ambient effects are neglected 
and the radial component of the stress tensor at the 
periphery is balanced by the pressure due to surface tension 
forces. However, the position of the jet surface is t iin- 

N. 

known. Since the radius-aeo is much larger than the perturb- 

s 

ation quantity- a, the ''"-boundary conditions are evaluated at 
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r . The location of this surface is given by the 
relation 




(5.-5,). 


where the operator- D_ is the usual substantial 

Dt 

derivative. Thus, the boundary conditions are expressed as 


B.C.l — 0 X (5.6) 

B.C.2 [~-p + "C^r] ( = (5.7) 


where p G , 



the surface tension pressure is given by 


a 



(5.8) 


B.C. 3 and 4 All velocities are finite at r=0 . 

s 

The rheological properties of viscoelastic materials 
m the linear region are often expressed m operational 
form as 


Pt = Q A (5.9) 

13 13 

where P and Q are polynomials with constant coeficients 

in some convected derivative operator such as the Olydroyd 

(65) or Jaumann (67) derivative. Linearization of the e 

convected derivative operator under steady flow conditions 

leads to polynomials in the operator d__ . P and Q may 

dZ 


then be written as 
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P = I + Woo • + A 2 V<J +• A m \C <3 l1_ (5.10) 

Q = 7o (‘ + /AiN / «>.|l-+/aiV« |L + ' ' ' • + (5 ' 1;l) 

The viscoelastic fluid has zero shear viscosity- T^ 0 and 
the X y A 2 ' • • • • f Am” /^l' >^2' ••••'/An are properties 

of the material. 

All variables are assumed to be expressible m the 

form 


$(^ 2 ) = FCO^ ( 5 - 12 ) 

By equation (5.12), the constitutive equation becomes 


r H= 


(5. IB) 


y[(h) - 70 


1- 

(■*^ + 1 

*(-!)" 

(. a*, )^m( 5 .13a) 

i-( 

^)a, + a v r x>- — -nfi 

< ftooj 

i N 


The equations of continuity and motion are similarly modi- 
fied to 


L _ 

r 5>r (Aco 


— ^V”Jl y _ — ^ -+- 

Qoo dr 


(5.14) 




(5.15) 
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= &;? + ^(h)[i| F (r|S) + (^V i ] = 0 (5.16) 

Equations (5.15) and (5.16) are crossdif f erentiated to 
eliminate the pressure term. A stream function 'p is 
defined as 


= = 
r a? 


-ili V* = -LSL 


0» r 


r ar 


(5.17) 


With the introduction of the stream function, a partial 

differential equation m 'p is obtained with the complex 

Reynolds number- Reoc (h) - 2acoVoep as parameter. 

^(h) 

D[ D + h R e oP (h) ] y = 0 
a oo 2 

where D = jh - jt ^ 4 f*. 

The complete solution to (5.18) is 

? = [V ! (r) +y 2 (r)]e~ hz / a “ (5.19) 

y?here y (r) is the solution to D'P 1 (r)=0 and 

y 2 (r) is the solution to (D + h Re^th) ] y , (r) =0 

2 
a oa 

The result is 


(5.18) 

(5.18a) 


y 


[ Ahj^_ /hr/aoo) + BhY^ (hr/a«, ) + CmJ-^ (mr/aj + (5.20) 
+ EmYj_ (mr/a«) )]e“^ z / aes 
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where m 2 =h 2 +h*Re DO (h)/2 (5.21) 

J n ( ) and Y n ( ) are modified Bessel functions of order n. 

The constants A,B,C,E are to be evaluated by the boundary 
conditions . 

Insertion of the boundary conditions, after appropriate 
modification with the stream function, results m a charact- 
eristic equation for the damping coefficient -h as a function 

i 

of the complex Reynolds number and a surface tension para- 
meter-S (h) = s' a«> /-»^>\) . The result is 



~ mJ Q (m) 1 + 2(m 2 -h 2 )h - (m 2 +h 2 )J 0 (h) =S(l+h 2 )h (5.22) 

. J i(m) J Ji(h) 


5.2 SOLUTION OF THE CHARACTERISTIC EQUATION 

Equation (5.22) permits an infinite number of solutions 
and any perturbation quantity wauld be given by the sum of 
an infinite series. At large axial distances, the profile 
rearrangement should be controlled by the lowest eigenvalue 
since its decay rate is the slowest. 

The lowest eigenvalue for the profile decay of a 
capillary jet of a Maxwell model fluid has been calculated 
by trial and error procedure and is given m figure 5-2 as 
a function of the final Reynolds number , with the Deborah 
number as parameter. For this viscoelastic fluid, the 
following definitions apply. The viscosity -V 7 (h) is given by 





' Y [(>) = 


(5.23) 




Figure 5-2 The lowest eigenvalue for a relaxi ng Maxwell model jet as cx function 
of the final Reynolds number, surface tension parameter 8-0, 

complex viscosity - / ^Ch') = J jj - h] * 


ij? 
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The final Reynolds number and Deborah number are expressed 
as 

(Re-co) = (Deb co ) = AO& L. (5.23a) 

(W) u*o 

Over the range of Reynolds numbers investigated, variations 
in the surface tension parameter between S=0 and S=60 did 
not affect the eigenvalues shown in the figure to any 
significant extent. 

The results of figure 5-2 show that a Maxwell model 
fluid relaxes faster than a Newtonian fluid of the same 
zero shear viscosity. This prediction is compared to 
experimental data in the next section. Under typical 

laboratory conditions of Voq = 200 cm/sec and (2Q*>) = 

-4 

0.2 cm., a relaxation time of 10 seconds is required 
for a Deborah number of 0.2 . This value is in agree- 

ment with estimates of the relaxation time of polymer 
solutions at high shear rates as calculated from independ- 
ent experimental data (55) . 

5.3 EXPERIMENTAL PROCEDURES 

Photographs were taken of capillary gets of Newtonian, 
non-Newtonian inelastic and viscoelastic fluids along their 
length in the region immediately outside of the capillary. 
Measurements were then made of the jet diameter as a 
function of axial distance which will be related to the 
rate of decay of stresses generated during the previous 
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shear flow. 

Photographic magnifications of approximately 10/1 
were obtained on a 4"x5" bellows camera equipped with a 
Polaroid back. A reference rod of known diameter which 
was sized to fit snuggly within the nozzle bore was also 
photographed at each location. This provided absolute 
dimensions to correct for variations m the photographic 
magnification along the jet length due to the unavoidable 
non-parallelism between the plane of the jet and the plane 
of the camera. A profile projector equipped with a micro- 
meter stage measured jets diameters from the photographs 
and the values reported are considered accurate to three 
significant figures. Chapter 3 provides a more detailed 
description of the procedures and equipment used . 

5.4 EXPERIMENTAL MEASUREMENTS OF PROFILE DECAY 

Measurements of the jet diameter as a function of 
axial distance from the nozzle exit were obtained for 
Newtonian, viscomelastic and viscoelastic liquids. A 
single nozzle diameter was used and data were taken over 
a range of fluid velocities. 

Newtonian fluids were included in order to verify the 
procedures and measuring techniques used m this investigat- 
ion with published results. For this purpose, mixtures of 
glycerin and water were employed as test fluids. Their 
properties are given m standard handbooks and have been 
summarized in figure 4-1. 
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The shapes of capillary jets of glycerin-water 
mixtures are shown m figure 5-3 as a function of the 
axial distance from the nozzle for a series of Reynolds 
numbers. All the jets are seen to undergo contraction 
within the range of variables explored . The length 
required for complete profile relaxation is seen to 
increase with the Reynolds number. 

The diameter ratio- ~)L = a M /a 0 , is plotted m 
figure 5-4 and is seen to be a decreasing function of 
the Reynolds number. Only at the highest Reynolds number 
does 'X approach the theoretical limit (35) of 'X = \f3/2 . 
Data points from this investigation fall smoothly within 
the dotted lines representing the scatter m published 
data (26,32,60). 

The prediction of exponential decay at large axial 
distances was tested by plotting In ( a ,) vs. In 
(z/a oo) as m figure 5-5 and the resulting straight lines 

confirm this expectation. From the slope of the lines, the 
damping coefficient is obtained. These values are plotted 
in figure 5-6 as a function of the Reynolds number. Also 
included m this figure are the experimental results of 
Goren and Wronski, indicated by the shaded area, and the 
theoretical solution for the lowest eigenvalue with the 
surface tension parameter- S=0. The profile decay rates 
found m this research are somewhat higher than those 
reported by Goren and Wronski. 
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AXIAL DISTANCE - 2/a 0 


Pigare 5-3 Diameter of capillary jets of glycerin- 
water as a fan ction of the axial 
distance from the nozzle and the 
Reynolds number^ 
nozzle diameter- 0.228 cm. 


DIAMETER RATlO-%^ 
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REYNOLDS NUMBER- Re 0 = / 2a 0 V 0 /\) 

» 

Figure 5-4 Diameter ratio of Newtonian jets 

as a -function ot the Reynolds number^ 
— — scatter in published data (26,32, 60), 
f!uid*-< 3 iycerin-wafer mixtures, 
hozxle diameter-0,228 cm, 



h - 2.7 


RSoo - A7. 1 


h =0.85 


AXIAL DISTANCE - Z/ao 

EjgureS-B Exponential profile decay in capillary jets of glycerin-water for a. 
series of Reynolds numbers} nozzle diameter- 0.228 cm. 
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AXIAL 01 STANCE- ^/a 0 
Fj0ure5-5 continued. 
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Figure 5-6 Damping coefficient for Wewfonian jets 
as a -function of the final Reynolds 

number; — Perturbation theory for 

S = 0 ? shaded area- data of Goren and 

WronsKi (32} , data of the current 

investigation. 
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The generally good agreement between the current 
measurements and published data, as typified by the 
comparisons m figures 5-3 through 5-6, are connotative 
of satisfactory techniques and procedures. 

The behavior of a vis coinelastic fluid undergoing 
stress relaxation was investigated by obtaining profile 
measurements for a relaxing capillary jet of 0.6% Carbopol. 
Figure 5-7 shows the shape of the jet as a function of the 
axial distance from the nozzle for three different Reynolds 

r • 

numbers-Re 0 , where Re 0 is the generalized non-Newtonian 
Reynolds number defined by Metzner and Park (53) as 

n 2-n 

Re = ( 2a d> < v o> © (5.24) 

Kgc ^ 1 

The constants K and n are the parameters m the power law 
constitutive equation (1.14) 

The inclusion of the Reynolds number as the primary 
parameter influencing the profile decay rate in viscom- 
elastic jets is based on its similar significance for 
Newtonian jets. However, some commentary is necessary 
regarding the particular choice of Reynolds number. If 
it is assumed that the 0.6% Carbopol regains its structure 
instantly, then the Reynolds number based on its zero shear 
viscosity would be the correct correlating group. With such 
a choice, the calculated Reynolds numbers for all three 
flow rates would be less than unity and it is unlikely that 
the behavioral differences evidenced in figure 5-7 would 
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AXIAL DISTANCE - 1 /a 0 

Figure 5-7 Diameter of cap il(av-y jets of 0.6%Carbofoo( as a function of the cocial 

distance from the nozzle and the Reynolds number:, nozzle diameter =0,7.75 cm 
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occur within the narrow Reynolds number range of 0<Re Q <l^ 

I 

But the use of the generalized Reynolds number-Re Q arranges 
the data into results that are similar to those published 
for Newtonian fluids m that expansion takes place for 

I ( 

Re =12.8 while contraction is" obtained for Re =51.0 and 
o o 

i « 

for Re o =90.2„ . Since Re Q is essentially a Reynolds number 
evaluated at the shear rate of capillary flow, , its usage 
implies that 0.6 % Carbopol does not exhibit any significant 
viscosity buildup above that possessed at the capillary 
exit. This could occur if the structural formation time 
was large compared to the time required for stress relax- 
ation. 

The jet diameter of 0.6 % Carbopol at a Reynolds 

I 

number of Re 0 =12.8 initially increases and then decreases 
before finally achieving a constant value with increasing 

t i 

axial distance. At Reynolds numbers of Re o =51.0 and Re Q = 
90.2 an opposite behavior is observed. The latter response 
has been reported in Newtonian jets (32) at Reynolds numbers 
between Re Q = 12-17 . 

The theory presented earlier in the chapter had pre- 
dicted that a capillary jet of a linear viscoelastic fluid 
would relax faster than a Newtonian jet of the same zero 
shear viscosity. In attempting to verify this prediction 
experimentally, it is important that the test fluid satisfy 
the assumption of a linear viscoelastic fluid. Practically 
speaking, this requires the polymeric, solutioft.-be '.sufficient 
ly dilute so that the ratio of (normal stresses/shearing 
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stress) is low. However, it has been found that profile 
decay measurements cannot be performed satisfactorily with 
dilute, low viscosity fluids due to the rapid growth of 
disturbances on their surface. 

The fluid chosen- 0.7% SCMC represents a compromise 
between these two requirements. Metzner et.al. (52) have 
reported that this fluid exhibits no measureable normal 
stresses at the shear rates encountered m a rheogonio- 
meter although normal stress effects are evidenced at the 
higher shear rates of capillary flow. 

The shapes of capillary jets of 0.7% SCMC are shown 

m figure 5-8 as a function of the axial distance from the 

» 

nozzle exit and the Reynolds number. Next to each curve 
is the Reynolds number based on the zero shear viscosity 
of 0.7% SCMC and the generalized non-Newtonian Reynolds 
number evaluated at the capillary flow conditions . Only 
at the lowest Reynolds number did the jet diameter become 
constant within the axial distances measured. Contraction 
also occured at Re 0 =13.6(with ReQ=305) although these data 
are not shown due to gravity effects at the low flow veloc- 
ity. 

Comparison of this data with the Newtonian results 
presented m figure 5-3 shows that a longer axial distance 
is required to reach a constant column diameter for a 0.7 % 
SCMC jet than for a Newtonian jet of the same zero shear 
viscosity. This is in direcr contradiction to the predict- 
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Figure 5*8 Diameterof capillary jets of 0.7% SCMC as a function of 

the axial distance from the nozzle and. the Reynolds number; 
nozzle diameter- 0.778 cm 


no 
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ion from linear theory presented earlier. The comparison 
could not be made m terms of the damping coeffiaent-h as 
values for the final relaxed diameter were not available 
for the 0.7% SCMC jets. 

If it is assumed that the structural formation time of 
0.7% SCMC is large compared to the time required for stress 
relaxation, then the decay rate should be a function of 

I 

some shear dependent Reynolds number such as Re Q . At 

t 

the values of Re Q shown m figure 5-8, it is unlikely 
that this mechanism is followed since the jet surface 
would show the effects of air friction and turbulence. 

The behavior evidenced by the 0.7 % SCMC jets would 
be consistent with a structural formation time which is 

t 

intermediate between the two cases considered of 1. an 
* instantaneous 1 regeneration of structure and 2 . an 
'infinite' time for viscosity buildup. 

An alternate explanation, originally proposed by 
Gill and Gavis is that a slowly decaying axial tension 
exists along the jet. Continuity arguements show? that 
this would cause the jet to expand. If such an effect 
is superimposed upon a jet contracting due to profile 
rearrangment, the net result would be to reduce the rate of 
profile decay. This explanation is m accord with the 
observed appearance of a viscoelastic jet undergoing breakup, > 
which is discussed m chapter 6. It will be shown that axial 
tensile stresses are present that cause large axial deform- 
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ations . 

5.5 SUMMARY OP EXPERIMENTAL RESULTS 

Measurements were obtained for the shape of capillary 
jets of glycerin-water mixtures (Newtonian) , 0.6 % Carbopol 
(viscomelastic) and 0.7 % SCMC (viscoelastic) . The Newton- 
ian results compared reasonably well with published data. 
Capillary jets of 0.6 % Carbopol expanded at low flow 
velocities and contracted at higher flow velocities. In 
order to explain the magnitude of this behavior, it was 

i 

necessary to assume that the structural formation time of 
0‘. 6 % Carbopol was large compared to the time required for 
stress relaxation. Inelastic fluids may therefore ^possess 
a relaxation time m the manner of viscoelastic fluids and 
the effects of this property can be exhibited m an exper- 
iment where the relaxation time and the deformation time 
■ 

are of comparable magnitude. 

Capillary jets of 0.7 % SCMC were found to decay slower 
than a Newtonian jet of the same zero shear viscosity. One 
possible explanation, which is m accord with the experiment- 
al data of chapter 6, is that a slowly decaying axial stress 

\ 

exists that serves to reduce the rate of profile decay. 

The profile decay data for both the viscomelastic and 
the viscoelastic jets have demonstrated that the relaxation 
of a long chain molecule is not necessarily m equilibrium 
with the relaxation of dynamic stresses. 



113 


CHAPTER 6 

BREAKUP OF A LAMINAR CAPILLARY JET 
OF A NON-NEWTONIAN LIQUID 

This chapter presents the experimental results 
obtained for the laminar breakup of capillary jets of 
Newtonian, non-Newtonian inelastic and viscoelastic 
liquids. A theoretical model developed by my colleague. 
Dr. J. Yerushalmi, for the stability of laminar visco- 
elastic jets is included and the use of his analysis 
is gratefully acknowledged. 

6.1 SOLUTION OF THE DYNAMIC EQUATIONS 

Consider an incompressible liquid jet emerging 
horizontally from a capillary and undergoing complete 
relaxation to a radius-ac© m a distance which may be 
neglected m comparison to the breakup length. The 
initially relaxed column then moves at a constant 
velocity-^, and is unaffected by gravity or the presence 
of the surroundings. Due to the presence of external 
vibrations or microscopic burrs within the nozzle, an 
infinitesimal disturbance is superimposed upon the motion 
of the jet and propagates as an axisymmetric wave. Thus, 
the jet is circular in crossection at all times and either 
expands or contracts. 

A cylindrical coordinate system is chosen (r, 0 , z) 
which moves with a constant velocity-V,*> „ Since the jet 
is taken to be initially relaxed, all components of the 
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velocity, stress tensor and pressure m the equations of 
motion and continuity represent small perturbation quanti- 
ties only. Non-linear terms are therefore second order m 
magnitude and are omitted. In the usual notation, the 
linearized equations of continuity and motion are given by 


i_ BCrV r) , etV& _ 
r Br dZ 


O ? jVr 

= - ££ + 

1 B CrTr< ) 

Cea j_ 

K at 

Br 

r dr 

r 32 



1 oCrT>e) 

. STTz-g; 

? St 

- + 
■dZ + 

r B2 

3 2 


The jet radius at any point may be expressed as 


( 6 . 1 } 

(6.2a) 

(6.2b) 


r = Q«, + $ (6.3) 

where (z,t) is the infinitesimal displacement of the 
free surface from the unperturbed radius. 

At the jet surface, the shear stress- T7 rz is zero 
since ambient drag effects are neglected and the radial 
stress tensor is balanced by the pressure due to surface 
tension forces. The difficulty m applying the boundary 
conditions is that the position of the jet surface is 
unknown. If the displacement ^<^3^ , the boundary con- 

ditions may be approximated as 

B.C. 1 (Trj) - O (6.4) 

r^eKso 


B.C. 2 (-pi- Trr ) = 


(6.5) 
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•where the surface tension pressure-p ff is given by 


Rr 





( 6 . 6 ) 


and S' is the coefficient of surface tension. 


B.C. 3 and 4 All velocities are finite at r=0 
Biot [2) and Bland (4) have shown that the 
rheological properties of any incompressible viscoelastic 
material m the region of very small deformations can be 
described by the following constitutive equation. 


= e ^ + X e 






(6.7) 


where ^ and A are the components of the strain and 
rate of deformation tensors, respectively, and E, 7^>»,C r and 
^ ^ are all real and positive. The coefficients E and ~y]co 
account for the instantaneous elasticity and the long term 
viscous flow respectively. The summation represents a 
discrete spectrum of relaxation times- l/<^ r 

Application of the Laplace time transform to equat- 
ions (6.1) and <(6.2a), and (6.7) results m 


± a Cry r) . = 


r dv 




0 


( 6 . 8 ) 


+ *!«[§? 


3 2* J 


(6.9a) 
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A 

£<* Vg =■ 




A 



9 z Vg ~| 


(6.9b) 


where <X is the transform variable and of (<x.) 
complex viscosity of eq„ (6.7) . The complex viscosity 
represents the linear response of any viscoelastic 
meterial and is given by 




+ 


Z 

i 




( 6 . 10 ) 


The boundary conditions are similarly transformed to 


B.C. !■ (,§f 


3r ) 


O 


B.C. 2' (-p 

/N 

c* ^ 


vis a. oo 

\ c? Vr~ 


(6.11) 

(0,K - <«• 13> 


'rs- A.OO 


Hereafter, the solution is identical to the method of Weber 
with the constant Newtonian viscosity replaced by the 

A 

complex viscosity . 

A 

A stream function ¥ is defined as 


c = -L ~ 

v r r 3? ^ v 2 


_L^i 
r dc 


(6.14) 


If equations (6.9a) and (6.9b) are crossdif ferentiated to 
eliminate the pressure and the stream function is introduced. 
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a linear differential equation in ^ is obtained as 
[d-^-]D?*0 (6.15) 

where D is a linear operator given as 


D = 


_ I b . 

arz. r 3r 



(6.16) 


A 

The stream function S* is assumed to be of the form 


?(r,Z) - [$,(0 + ^(r) ] . {6 * 17) 

where k, the wave number, is related to the wavelength 
by k =2 IT*/ £ . 

is the solution to D 'f ' =0 
and 4 / 2 is the solution to [ D — <x q/y\ fc)] = o 

A 

The complete formal solution for is given by 

y(r,z)= r jc 1 I 1 (kr)+C 2 I 1 (gr)-K: 3 K ;L (kr)“fC 4 K ;L (gr)J x (6.18) 
x 

where g I 2 = k 2 + °< ^ /-yj (*) 

I ( ) and K n ( ) are modified Bessel functions of 
order n and the 0 ^, 02,03 C 4 are constants to be determined 
from the boundary conditions. 

_ Application of the boundary conditions leads to the 
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characteristic equation relating o<. and k 

« 2 + 2-y (-Ok 2 [ I,' (ka«,) -2kg_ I, (ka«,) (ga„) = (6.19) 

« I c (ka„) - Il<ga „) J 

= <s"k (l-k^a.?) I^(kaco) a 2 -k 2 

where I^( ) denotes differentiation of I n ( ) with respect 
to the arguement. 

The solution to the characteristic equation determines 
the stability of the capillary jet. Depending upon whether 
the real part of <x. is positive, zero or negative, the 
initial disturbances will be amplified, neutrally stable, 
or dampened, respectively. 

6.2 SOLUTION OF THE CHARACTERISTIC EQUATION 

A considerable simplification of eq. (6.19) results 

when the wavelength is much greater than the jet radius, 

O-po . or, equivalently, when (kao©) 1 and when((gaoo) 

1. This condition is fufilled by low speed capillary 

jets of Newtonian liquids. Under these circumstances, the 

Bessel functions may be approximated by the first term of 

their series expansions and can be written as 

I (ha*,) 1 ,* l 1 (ka 00 )=I 1 (ga 0 o) 

° ( 6 . 20 ) 

I-^ka*,) kaco/2 ; x-l (ga^) 'ga^ /2 

With these simplifications, the characteristic equation 
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becomes 


S^d-K^S.) =o 


The capillary jet will be destroyed by the fastest growing 

* A* 

wave- <x which is obtained from the solution of =0 

and is given by 


&C * - 


r 2. oak 


2. + 3*jfr0 . 1 2. p 

^ V 


( 6 . 22 ) 


The highest growth rate of any fluid will be for an mviscid 

y\ ^ 

liquid in which 'M (°0 =0 . Corresponding to Cx is the 
* 

wave number-k 




(6.23) 


The stability of a capillary jet will now be 
considered m relation to the properties of the complex 
viscosity. For the convenience of the reader, the general 
form of ‘V'jW is repeated. 




*+* ^oO 1* 


- 

Z C r 
b r t 


( 6 . 10 ) 


CASE 1 NEWTONIAN LIQUID 


For _ a Newtonian liquid of constant viscosity - 7 ^ , 

equations (6 .22) and (6.23) become 
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Z flao 


■4* Z 0.t» 



(6.24) 


(6.25) 


* * 

where the subscript zero in <X and k Q denotes 
Newtonian values. Equations (6.24) and (6.25) are iden- 
tical to the results of Weber. 

The dimensionless wavelength for a capillary jet of a 
Newtonian liquid is given by 

If the breakup time is assumed to be inversely proportional 
to the growth rate, then the breakup length of a Newtonian 
jet can be written as 

L » CV^/oc* < 6 - 27 > 

where L is the breakup length and C is a constant of 
proportionality which must be determined experimentally. 

By virtue of equation (6.24) the dimensionless breakup 
length is obtained as 

Vn/ £, Wg.oa 

Re 


i + 


3U 

fz ecui 




(6.26) 


L 

2A«s 


c 


(6.28) 
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where We W; the Weber number is given by We.„=2 £ a ® Vw/c* , 

t 

and Re (?; the Reynolds number is defined as Rc.oo=2a&» ^ /y^ 0 . 

CASE 2 VISCOELASTIC FLUID WITHOUT STRUCTURE 

The relative stability of a capillary jet of a visco- 
elastic fluid can be compared to that of a Newtonian 
liquid by dividing eq. (6.22) with eq. (6.24) to give 


Z± 

2 + 




(6.29) 


* * 

The ratio (°<- 0 / ) thus depends upon the relative magni- 
tude of with respect to Since the 

complex viscosity of a viscoelastic fluid with E=0 is 

bounded by C°°) < '‘l 0**) < "‘j (<>) it follows from eq. 

* * 

(6.29) that c <. o /c« <1 and, therefore, a capillary jet of 

such a fluid is less stable than a Newtonian jet of the same 

zero shear viscosity- ' v \(o) -^o . 

2 

A plot of c< vs. k , derived from eq. (6.21) for 

A. . 

a viscoelastic fluid with finite yields a curve which 

lies above the curve corresponding to a Newtonian fluid of 
the same zero shear viscosity and below the curve corres- 
ponding to an inviscid fluid =0. All three 

2 

cross the k axis at the same points. 


curves 
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CASE 3 ELASTIC SOLID 

When E is finite with and all the r equal to 

zero, the material behaves as an elastic solid for which 


Of 

It can be shown that if 



(6.30) 


E c (l - K x dw) 

( o C\ oo 


(6.31) 


* 

there will always be a real positive eigenvalue of <X 
causing instability while if inequality (6.31) is reversed, 
the jet is unconditionally stable. 


CASE 4 VISCOELASTIC FLUID WITH STRUCTURE (E^O) 

The same criteria governing the stability of a jet 
of an elastic solid, case 3, can also be shown to be 
applicable to the case of a viscoelastic fluid with 
structure . 

The interested reader may consult the thesis of Dr. 
Yerushalmi (94) or the paper by Goldin et.aL. (29) for 
details of this development. 


6.3 DISCUSSION OF THEORETICAL RESULTS 

In summary, it has been shown that a laminar capillary 
jet of a viscoelastic liquid whose complex viscosity is 
bounded by ^ (o) *>1 0 , is more unstable that a Newtonian jet 
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of the same zero shear viscosity- . If, however, Tf(p<) 

is unbounded the jet may become unconditionally stable. 

It is worthwhile to discuss some of the assumptions 
leading to these conclusions. Firstly, the liquid column 
was taken to be completely relaxed. As noted m the 
previous chapter, the ]et , upon emerging from the capillary, 
undergoes a relaxation m which the stresses and velocity 
profile (of the capillary shear flow) decay with increas- 
ing axial distance. The rearrangement results in a relaxed 
jet diameter that is different from the capillary bore. 

For laminar Newtonian jets, it has been shown (32,57) by 
measuring the axial distance required for the jet to 
reach constant diameter, that the length required for 
complete relaxation may be neglected m comparison to the 
breakup length. Middleman and Gavis (61,62) and Metzner 
et.al. ( 52 ) have obtained similar data for jets of visco- 
elastic liquids but only report the ratio of the final 
relaxed diameter to the nozzle diameter. It may be 
inferred from their lack of comment as to the relaxation 
length that this distance was not unusual and was of the 
same order as for Newtonian jets. However, Gill and Gavis 
(27) interpret their experimental results to mean that a 
slowly decaying tensile stress exists m jets of visco- 
elastic liquids even at considerable axial distances. 

For viscoelastic fluids consisting of dissolved poly- 
mers, there is also a structural relaxation time for the 
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long chain molecule to achieve an equilibrium, spatial 
configuration with the locally changing shear rate. If 

the molecular relaxation time is of the same order as the 

-2 -3 

jet breakup time, about 10 to 10 seconds, the jet may 
never be fully relaxed . 

The linear stability analysis just described had 
assumed that the propagating disturbance was infini- 
tesimally small. On that basis, terms which were second 
order m magnitude were discarded from both the equations 
of motion and the rheological constitutive equation of the 
fluid. There is justification for neglecting the inertial 
terms m the equations of motion as the same assumption in 
the treatment of Newtonian jet breakup has produced theo- 
retical predictions m good agreement with experimental 
data. However, the use of a linear constitutive equation 
for viscoelastic fluids which exhibit measureable non- 
Imearities such as shear dependent viscosities and normal 
stresses is a much more restrictive approximation. The 
analysis may be applicable to very dilute viscoelastic 
fluids although Oliver (68) has measured normal stresses 
at concentrations down to the 0.01% level m jet thrust 
experiments. In the end, experimental data must be used to 
access the validity of the assumptions. As the photographs 
m the experimental section will show, marked non-linear 
phenomena are present at large axial distances even for 
jets of very dilute viscoelastic liquids. 
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6.4 EXPERIMENTAL RESULTS - QUALITATIVE OBSERVATIONS 

Marked differences were found to exist between the 
behavior of capillary jets of Newtonian, non-Newtonian 
inelastic and viscoelastic fluids. The distinctions are 
summarized by the photographs of figure 6-1, which shows 
typical appearances of the jets formed by each class of 
fluids . 

In laminar jets of Newtonian liquids, an infinitesimal 
disturbance generated within the capillary is propagated as 
an exponentially growing wave with a constant wavelength, 
according to equations (6.24) and (6.26). Capillary jets 
of low viscosity liquids such as 75% glycerin-water and 
ethylene glycol, which are seen m photographs (a) and (b) 
respectively, clearly illustrate this mode of instability. 
Viscous Newtonian jets are also broken up by the regular 
growth of a disturbance wave but, as shown m photograph 
(c) for a 97% glycerin-water jet, the wave character is 
not as uniform. Just before breakup, these high viscosity 
jets may also form fine threads over relatively short 
lengths . 

The non-Newtonian, inelastic test fluids possessed 
certain rheological similarities. They all displayed a 
yield stress and a shear dependent viscosity but did not 
exhibit any of the phenomena normally associated with 
viscoelasticity, e.g. a Weissenberg effect. 

Laminar capillary jets of 0.1% Carbopol and the MPA 
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60 dispersions, photographs (d), (e) and (f) respectively, 

appear to have characteristics similar to those of Newtonian 
jets m that they break up due to the growth of a disturb- 
ance wave. When the Carbopol concentration has increased 
to 0.2% from 0.1%, the fluid becomes considerably more 
viscous and the capillary jet formed by 0.2% Carbopol, 
shown m photograph (g) , shows a barely perceptible wave 
propagation. 

A 0.6% Carbopol solution has the consistency of a 
semirigid gell and the capillary jet of this fluid, seen 
m photograph (h) , is unstable but only at very high 
velocities. No wave characteristics are apparent on the 
liquid column. The drag force of the surrounding air is 
no longer negligible under these high velocity flow condi- 
tions and imparts an irregular motion to the capillary jet. 

An aqueous dispersion of silica at a 12.4% concentra- 
tion is an inelastic liquid with a high yield stress. It 
forms capillary jets that show no regular wave pattern. At 
the intermediate velocity at which photograph (l) was taken, 
the column is stable and the “screw-like" disturbance 
pattern that is seen on the liquid column does not grow. 

The appearance of the capillary jet is remindful of melt 
fracture in polymer extrudates. When high velocities are 
reached, the jet becomes turbulent and breaks up rapidly 
due to the growth of non-axisymmetnc disturbances, as 
shown m photograph (j) . The high velocity silica jet 
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resembles the turbulent breakup of a low viscosity Newtonian 
liquid despite the fact that, at rest, the 12.4% silica 
dispersion is a stiff gel. 

The viscoelastic test liquids used xn these experi- 
ments all possess shear dependent viscosities, normal 
stresses (68,77) and show marked drag reduction at very 
low levels of polymer (23,36,76,84). At the same con- 
centration, Separan and Polyox solutions display stronger 
viscoelastic effects than those of SCMC. 

Dilute viscoelastic solutions of SCMC m water form 
jets which initially display a growing wave with a clearly 
defined wavelength, as m photograph (k) . The growth of 
the wave is, however, arrested before breakup occurs and a 
series of droplets connected by threads appears. These 
•Threads thin with distance and eventually lead to the breakup 
of the liquid column. 

In moderately elastic jets such as 0.25% Separan in 
water, shown in photograph (1), no wave motion is discern- 
ible and the first visible disturbance appears as a large 
droplet, isolated m space from any systematic growth 
pattern. At greater axial distances, the jet is composed 
of a series of random sized droplets connected by random 
lengths of threads, as illustrated by photograph (m) of 
0.25% Polyox. 

Dilute viscoelastic jets of 0.05% Separan, seen in 
photograph (n) , exhibit an intermediate behavior between 
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0.25% SCM.C and 0.25% Separan or 0.25% Polyox. The 
disturbance initially propagates as a wave of irregular 
amplitude with a wavelength that increases with distance. 
Gradually, the waves form a string of droplets connected 
by threads as m the case of the more elastic jets. The 
droplet distribution, while not completely random, is less 
regular than m 0.25% SCMC. 

Capillary jets of viscoelastic liquids are thus dis- 
tinguished by the formation of droplets connected by threads 
which thin with distance and ultimately break. Though 
threads have been observed near breakup for viscous 
Newtonian jets, they exist only over relatively short 
distances. The large axial lengths over which the droplet- 
thread configuration is propagated appears to be a character- 
istic peculiar to viscoelastic jets. 

Since both non-Newtonian inelastic and viscoelastic 
liquids have shear dependent viscosities, the considerable 
difference m the appearances of their respective jets 
undoubtedly is associated with the presence or lack of 
elastic properties. 

The complete breakup of a capillary jet of ethylene 
glycol is shown m the sequence of photographs of figure 
6-2. The distances quoted in the legend represent the 
length from the nozzle tip to the midpoint of the photo- 
graph. Photograph (c) of this figure clearly shows the 
wave character of the propagating disturbance and also 
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depicts the well defined wavelength. The uniform distance 
between free droplets is an additional indication of the “ 
wave character of the propagating disturbance. 

Figure 6-3 shows the breakup sequence for a jet of 
97% glycerin-water and is illustrative of the behavior of 
a high viscosity, Newtonian liquid. As shown m photograph 
(c), the disturbance initially propagates m the form of a 
growing wave. However, near breakup, photograph (d) shows 
how the uniformity of the wave is altered. Fine liquid 
threads often appear m this region but only over rela- 
tively short distances. 

The wave propagation of disturbances seen on capillary 
jets of 0.1% Carbopol and of 24% MPA 60 in xylene, figures 
6-4 and 6-5 respectively, is strikingly similar to the 
appearance of low viscosity Newtonian jets. A related 
behavior is shown by 0.2% Carbopol jets in figure 6-6 even 
though the amplitude of the disturbance wave is consider- 
ably diminished. The Carbopol solutions and the MPA 60 
dispersions exhibit definite yield stresses and show high 
apparent viscosities at low rates of shear. Yet all 
atomize readily and the breakup lengths are not indicative 
of any appreciable stability. 

Solutions of 0.6% Carbopol are semi-rigid gells when 
at rest with a moderate yield stress. Capillary jets formed 
by these solutions are unstable, but only at very high 
velocities. Figure 6-7 shows the complete breakup 
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sequence for a capillary jet of 0. <3$ Carbopol moving with 
an average velocity of 4710 cm/sec. Air resistance at 
these high velocities produces the non-axisymmetric motion 
of the liquid column that ik evident m photographs (c) 
and (d) . The wave propagation which characterized the 
capillary jets of 0.1% Carbopol, 0.2% Carbopol and the 
MPA 60 dispersions is absent with 0.6% Carbopol. This is 
also confirmed by the random nature of the particles 
formed after breakup, as shown m photograph (e) . 

The sequence of photographs appearing m figures 
6-8 through 6-10 shows the changing behavior of a 
capillary jet of a 12.4% silica dispersion. Photograph 
(a) of figure 6-8 depicts the formation of a non- 
axisymmetric disturbance on the liquid column as it emerges 
from the nozzle at an average velocity of 1080 cm/sec. 

The amplitude of the disturbance at a distance of 20.3 cm 
from the nozzle, photograph (b) , does not appear to have 
increased over that at the capillary exit and actually 
seems to have dampened at a distance of 35.6 cm., as m 
photograph Oc) . The jet was stable within the 5 foot length 
available for photography. When the jet velocity has 
reached 1950 cm/sec, the montage pictured m figure 6-9 
shows large amplitude disturbances visible on the surface 
of the jet which eventually result m its destruction. The 
growth of these disturbances appears to be random and this 
is further evidenced by the non-uniform nature of the 



particles formed at breakup, as shown in photograph (d) . 

At the highest jet velocity, 2850 cm/sec, the 12.4% 
silica jet, seen m the sequence of figure 6-10, undergoes 
turbulent breakup and is completely fragmented into indi- 
vidual particles within a relatively short distance. In 
none of this series of photographs is there evidence of 
the wave propagation which characterized the capillary 
jets of Newtonian liquids, 0.1% and 0.2% Carbopol solutions 
and the MPA 60 dispersions. 

It is surprising to observe that 12.4% silica jets 
atomize so readily especially when one considers that 
this inorganic dispersion is a stiff gel, when at rest, 
with a high yield stress. Anothermteresting observation 
is that the lowest velocity silica jet is stable whereas 
the higher velocity jets are unstable. These phenomena 
will be examined at greater length in the discussion of 
experimental results, section 6-6. 

Solutions of 0.25% SCMC are very slightly elastic and 
the jets formed by this fluid initially show a growing wave 
with a clearly defined wavelength. The growth of the wave 
ceases before breakup occurs and a string of droplets 
connected by thin threads is formed. The droplets are found 
at regular distances from one another indicating their 
development from a wave of constant wavelength. This 
behavior is indicated m the sequence of photographs of 
figure 6-11. Secondary instabilities in the form of very 
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small droplets may also be generated on the threads, as 

i 

demonstrated in the last photograph of this figure. 

In capillary jets of 0.25% Separan and 0.25% Polyox, 
which are considerably more elastic than 0.25% SCMC, 
wave characteristics are absent and the first visible 
disturbance appears as a large droplet, isolated m space 
from any systematic growth pattern. The thread lengths 
connecting the droplets are , randomly distributed and thin 
with distance, eventually leading to the breakup of the 

t 

column. A series of photographs taken along with the length 
of a 0.25% Separan jet illustrates these phenomena and is 
presented m figure 6-12. The behavior of a 0.25% Polyox 
jet is similar to that shown by 0.25% Separan. At high 
velocities, the threads are able to undergo large ampli- 
tude, three dimensional disturbances without breaking, as 
shown in photograph (a) of figure 6-14. At a concentrat- 
ion level of 0.5% Separan, no observable disturbances were 
detected within the limits of the 5 foot distance which 
could be conveniently photographed. 

The photographs of figure 6-13 show how a 0.05% 

Separan jet changes m appearance from an initial wave-like 
contour to a configuration of droplets connected by threads. 
The threads are seen to thin with increasing distance from 
the nozzle and, at the same time, the droplet diameter 
increases. Jet breakup results from disruption of a thin- 
ning thread. In this dilute viscoelastic jet, there is a 
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preponderance of 'twxn' droplets xn juxtapositxon and 
separated from the rest by long threads. Photograph (d) of 
thxs fxgure xllustrates the xnterestxng symmetry of the 
'twxn 1 droplets. 

Another type of instabilxty xn whxch the lxquxd column 
possesses a screw orxentatxon has been observed wxth dxlute 
vxscoelastxc jets. There xs only a narrow range where thxs 

i 

phenomena xs dxstxnctly vxsible. Photographs (b) and (c) 
of fxgure 6-14 shows a 0.05% Guar Gum jet travelxng at an 
average velocxty of 350 cm/see and formxng a stable helix 
whose pxtch xncreases wxth dxstance from the nozzle. Near 
breakup, a dxsturbance wave wxth an xrregular, but identi- 
fxable wavelength can be seen superxmposed upon the helxx. 
Another example of a screw instabilxty can be seen xn photo- 
graphs (d) through (x) for a 0.05% Separan jet moving with 
an average velocxty of 690 cm/sec. The dxsturbance forms 
close to the nozzle but later dampens out and fxnally 
reverts to the droplet- thread conf xguratxon whxch has 
characterxzed the vxscoelastxc jets studxed xn thxs 
xnvestxgatxon. The screw xnstabxlxties observed for the 
capxllary jets of 0.05% Guar Gum, 0.05% Separan and the 
12.4% aqueous sxlxca dxspersxon resemble the shapes assumed 
by molten polymers undergoing melt fracture. 

6.5 SUMMARY OF QUALITATIVE EXPERIMENTAL RESULTS 

A short summary wxll be present xn order to provxde a 
concxse descriptxon of the qualxtative phenomena just 
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described . 

In capillary jets of low viscosity Newtonian liquids, 
disturbances are propagated as a growing wave with a 
constant wavelength. The free droplets formed after break- 
up are uniformly spaced and are an additional indication 
of their creation from a wave pattern. High viscosity 
Newtonian jets behave similarly except that, at breakup, 
the regularity of the wave is altered. Fine threads of 
liquid often appear m thas region but only over relatively 
short distances. 

The less gellified, viscomelastic fluids formed 
capillary jets which resembled the appearance of Newtonian 
jets m that a disturbance was propagated as a growing wave 
with a constant wavelength. Capillary jets of the highly 
gellified, viscomelastic fluids were unstable only at high 
velocities and displayed no wave characteristics. The 
instability of all liquids m this category is surprising 
since, at the low shear rates associated with the growth 
of symmetric disturbances, the viscomelastic test fluids 
show high apparent, viscosities . Viscosity is a stabiliz- 
ing influence and yet the breakup lengths of capillary 
jets of viscomelastic liquids were not indicative of any 
exceptional stability. 

Capillary jets of viscoelastic liquids are character- 
ized by the formation of droplets connected by threads 
which thin with distance and ultimately break. Even for 
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very dilute viscoelastic jets, which, initially show a 
disturbance wave, the droplet-thread configuration is 
ultimately formed. Concentrated solutions of the more 
elastic polymers form capillary jets which show no wave 
characteristics. The first visible disturbance appears as 
a large droplet, isolated m space from any systematic 
growth pattern. 

The photographs m figures 6-1 through 6-14 have 
served to qualitatively demonstrate the differences 
between the capillary jets of Newtonian, non-NeWtonian 
inelastic and visco.elastic liquids. Section 6-6 is devoted 
to a more detailed and quantitative presentation of the data. 

6.6 QUANTITATIVE EXPERIMENTAL RESULTS 

The experimental data to be presented in this section 
involves measurements of the jet breakup length and the 
wavelength of the propagating disturbance wave, as obtained 
from still , photographs and motion pictures. Results for 
capillary jets of Newtonian, non-Newtonian inelastic and 
, viscoelastic liquids are given m subsections 6.6-1, 

6.6-2 and 6.6-3 respectively. 

Either pressurized nitrogen or the crosshead of an 

> 

Instron mechanical tester mov]ed the piston of a cylinder 

t 

containing the test fluid and forced the liquid through 
hypodermic needle tubing. Further details on the equipment 

i 

and procedures used to obtain jet breakup data are provided 
m chapter 3. Fluid properties are discussed in chapter 4. 
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6.6-1 CAPILLARY JETS OF NEWTONIAN LIQUIDS 

In an experimental investigation of hydrodynamic 
stability, it is important that the apparatus introduce a 
low level of disturbances into the liquid. Data were 
obtained with Newtonian fluids m order to ascertain 
whether the equipment and procedures used m this work 
produced results which were comparable to those of previous 
investigators. The test fluids consisted of water, ethylene 
glycol, 75% glycerin-water, 97% glycerin-water and pure 
glycerin. 

The dimensionless breakup length— L/2a 0 , based on the 
capillary diameter, was measured for the capillary jets 
formed by each of the test liquids and, m the laminar 
range, was found to increase with velocity and nozzle 
diameter. These results are plotted m figures 6-15 
through 6-19. Where the velocity range is sufficiently wide, 
as m the case of the water jet fl the L vs „ V curve goes 
throu^i a maximum as air resistance becomes significant. 

Further increases in velocity bring about a decreased 
breakup length. 

As noted in chapter 5, profile relaxation causes the 
diameter of the initially relaxed liquid column to change 
from that of the capillary bore. The experimental measure- 
ments were converted to values based on the relaxed diameter 
using published data (26,32,60) for ~)C = a„ /a as a function 
of the capillary Reynolds number. By continuity arguements. 
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Figure G-IS Breakup length of a. capillary jet of water a6 a. 

-function of fluid velocity' f loud fempera.fu.re- 
J ZB> c ’C } noxxte diameter- 0,02.63 cm. 
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Figure 6-l<3 Breakup length of acafai Mary jet of ethylene glycol 
as a function of fluid velocity and nozzle d iameter$ 
fluid teniperature-25°c J nozzjediamefef r : o-o.O'ZiG 
cm., A-0. 0263 cm, x-0,04f 4 cm. ; D-O.08G& cm. 
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AVERAGE FLUID velocity- V 0 cm/s®c. 


Figure G'(7 Breakup length ofacaphlary jet of 75%giycertn- 
waier a<5 a Function of fluid, vetoclty; fluid 
temperature- 25 °C , nozzle dia.rneter*Q.Ot c IO cm. 



AVERAGE FLUID VELOCITY' Vo cm/sec, 

Figure G-1& Breakup length of acapiilary jet of CJ7% glycerin' 
wafer ae> a. -function of fluid velocity, f lutd 
temperature -25 °C, nozzle diameter- 0.041 4 cm. 
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Figure <b‘\°> Breakup length of a capillary jet of 
glycerin 06 a function of fluid 
velocity 5 fluid, tern joe nature- 61 °C ? 
nozzle diameter- 0.0 1^0 cm. 
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the following relationships apply 

v =v 0 //^ 2 7 Re « = Re o/X ' We «r - We o / % 3 ( 6 - 32 > 

where a<x, , V OT , (Re),*, = 2aoo Voo /v and (We),^, = 

t 

2 

2. aoo / <S ate the relaxed jet diameter, velocity, 

Reynolds number and Weber number, respectively. The 

* / 

corresponding quantities relative to the capillary dia- 

j 

meter are designed by the subscript zero. 

The converted data were compared to the extensive 
results of Kroesser (39) and Grant (33), as correlated by 
Kroesser, for the breakup length of a laminar Newtonian 
jet. As shown m figure 6-20, the agreement is excellent, 
indicating that the equipment and techniques are 
satisfactory. 

The wave length of the propagating disturbance was 
measured directly from still photographs of the capillary 
jets near breakup over a range of nozzle diameters and 
fluid velocities. It is compared to the prediction of 
equation (6.26) m table 6-1. This relation, rewritten 
below for the convenience of the reader, showe that the 
wavelength will be a weak function of the velocity because 
of the dependence of the relaxed diameter- 23 .^ on the 
velocity. The good agreement obtained between experi- 
mental and theoretical values is an additional ^conf irm- 
aticn of the linear stability theory. 
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Figure G-*2.0 Correlation of Mewfonian jet breoKap 

data, o-wafer, A- ethylene glycol, 0-75% 
glycerin- water, n~°l7 0 /o glycerin- water, 

X-pure glycerin, Correlation of 

Kroesser and Middleman (33 ) 




Table 6-1. Wavelengths of Disturbance Waves in Laminar Capillary Jets of Newtonian Liquids 


Fluid 

Nozzle 

Diameter 

Average Fluid 
Velocity 

Experimental 

Wavelength 

Theoretical Wave- 
length-Eq. (6-26) 

Water 

0.0263 cm. 

190 cm. /sec. 

0.08-0.10 cm. 

0.10 cm. 



310 

0.10 

0.10 



522 

0.10 

0.10 

Ethylene Glycol 

0.0216 

283 

0.10-0.12 

0.11 



515 

0.10-0.12 

0.11 



888 

0.10-0.12 

0.11 


0.0263 

622 

0.11 

0.13 


0.0414 

318 

0.20 

0.20 



590 

0.18 

0.19 


0.0868 

250 

0.38 

0.39 



290 

0.35 

0.39 



590 

0.35 

0.38 



820 

0.35-0.40 

0.38 

75% Glycerin- 





Water 

0.0868 

270 

0.49 

0.45 

97% Glycerin- 

0.0414 

708 

0.47 

0.50 


Water 
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6.6-2 CAPILLARY JETS OP VISCOINELASTIC LIQUIDS 

The breakup length of 0.1% Carbopol is shown xn fxgure 
6-21 as a functxon of the average fluxd velocxty and the 
nozzle dxaraeter. As the jet passes from lamxnar to turbu- 
lent flow, the L vs. V curve exhxbxts the same character- 
xstxc shape as for a Newtonxan lxquxd . For lamxnar jets, 
the breakup length xncreases wxth both the average fluid 
velocity and the nozzle diameter. However, as air 
resistance becomes significant, the’ L vs. V curve goes 
through a maximum and further increases xn velocity result 
in a decreased breakup length. When the jet becomes 
turbulent, the breakup length again increases. 

The stability of 0.1% Carbopol was compared to that 
of water, which behaves as an xnviscid liquid under these 
conditions, by plotting the breakup length against the 
average fluid velocity for a 0.0263 cm nozzle. Figure 6-22 
shows that when both curves are linear, the stability of 
the ideal fluid is greater than 1 that of the 0.1% Carbopol. 
This result is surprising since the linear theory presented 
earlier predicts that viscosity is a stabilizing factor; 
yet the viscosity of 0.1% Carbopol is higher than that of 
water at all shear rates and their surface tensions are 
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AVERAGE FLUID VELOCITY-V 0 cm/sec. 



0 GOO 1200 1800 2400 3000 

AVERAGE FLUID VELOCITY- V 0 cm/sec 


Figure G-2| Breakup length of a. capillary jet of 0,i% 
Carbopol as a. funcf ion of f laid velocity 
and nozzle diameter; nozzle diameter: 
A'O.o^an., x- o.04i4cm. ; a-o.0868 cm. 
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Figure b-'L'l BreaKup length of capil lary jets of 

wafer and 0.\°/o Cartoopol as a function 
of fluid, velocity 5 o -water, x- 0. 1 °/o 
Carbopol, nozzle diameter- 0.02 63 cm. 
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roughly equivalent. A more dramatic comparison is shown m 
figure 6-23 where the breakup length of a laminar capillary 
jet of water is. seen to be 2% times greater than that for 
0.1 % Carbopol at the same average fluid velocity and nozzle 
diameter. However, as shown m figure 6-22, the viscous 
consistency of the 0.1 % Carbopol does inhibit breakup 
due to air resistance. The maximum m the L vs. V curve 
occurs at much higher velocities than for water. 

The dimensionless breakup length- L/2aoofor a capill- 
ary jet of a Newtonian liquid is predicted from equation 

h 

(6.28) to be proportional to (2a oo) for a low viscosity 

fluid and to be independent of (2a oo) for a high viscosity 
liquid . This result was confirmed by the experimental 
data of figure 6-20. In order to test the diameter 
dependence for 0.1 % Carbopol jets, the figure 6-21 data 
at an average velocity of 300 cm/sec, which is within the 
laminar range, were replotted as (L/2a Q ) vs. (2a 0 )^ in 
figure 6-24. The nozzle diameter, rather than the relaxed 
jet diameter, was used since no correlations for the 
diameter change are available for 0.1 % Carbopol. Based 
on the analagous Newtonian data, a maximum error of 13.4 % 
would be expected. The linear dependence shown m figure 
6-24 displays the mviscid nature of the 0.1% Carbopol jet. 

If the fluid jet is assumed to completely regain its 
zero shear viscosity after profile relaxation, the disturb- 
ance wave will grow at a rate corresponding to this 
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Figure G'24 Effect of nozzle diameter on 
loreaKup length of capillary 
jets of 0,1 % Carbopol ; 
average fluid velocity = $00 cm/sec. 
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viscosity since the shear rates of wave propagation are 
quite small. A solution of 0.1%, Carbopol shows moderate 
viscosities at low shear rates and yet the capillary jet 
of this fluid behaves as an inviscid liquid . This response 
is possible if the jet initially relaxes of all stresses 
originating from the previous capillary flow, but does not 
regain its zero shear viscosity because its structural 
relaxation time exceeds the jet breakup time. In that case, 
the behavior of the jet would correspond to that for a 

E 

liquid having a radially varying viscosity roughly equiva- 
lent to that existing at the exit shear rate. This 
prospect will be explored m the presentation of data for 
the other viscomelastic jets. 

The breakup length of capillary jets of 24% MPA 60 
in xylene and 24% MPA 60 m mineral spirits increases with 
the average fluid velocity and with the nozzle diameter, as 
shown m figure 6-25 and 6-26, respectively. All data 
points have been omitted for clarity because of the 
considerable scatter amongst them. Due to the thixotropic 
nature of the fluids, time dependent procedures such as the 
period spent m the hydraulic cylinder prior to a run may 
have an effect on the jet behavior. However, the tortuous 
path taken by the liquids through the fine pore filters, 
prior to entering the capillary, should have homogenized 
them. Another possibility for the variation is the presence 
of air bubbles introduced during the mixing that is required 
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Figure <2>*2G BreaKajD length of a capillary jet of '2.4% MPA GO in 
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to initially break the gel and pump the fluid into the 
hydraulic cylinder. 

There is no maximum observed m the L vs. V curves 
within the experimental range of velocities considered. 
Another difference from previously described phenomena is 
that the curves do not appear to pass through the origin 
as with Newtonian and 0.1% Carbopol data. This is 
consistent with the hypothesis that the fluid has not 
regained its zero shear viscosity after profile relaxa- 
tion but instead maintains a viscosity of the same 
magnitude and distribution as at the capillary exit. The 
explanation may be made clearer by recourse to the illus- 
tration of figure 6-27. The dotted lines represent the 
linear portion of the L vs. V qurves for capillary jets of 
three Newtonian liquids of viscosities < 7 ^ < y c , 
respectively. At velocity- V g/ the non-Newtonian inelastic 
fluid has an apparent viscosity and the breakup 

length- L a corresponds to that of a Newtonian jet with 
constant viscosity- . At some lower fluid velocity- Vj^, 

the average apparent viscosity of the viscoihelastic. fluid 
increases to T/m and the breakup length is 
approximately that for a Newtonian jet of constant viscosity- 
■y . If the argument ; is continued m like fashion, the 
extrapolated L vs. V line does not pass through the origin. 
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Figure 6-27. The behavior of a capillary jet 
of a viscomelastic liquid. 

A crossplot of the linear data m figures 6-25 and 6-26, 
taken at a constant velocity of 500 cxn/sec, is given m 
figure 6-28 and shows (L/2a 0 ) to be proportional to (2a 0 )% 
for both fluids. As with the 0.1% Carbopcl data, the nozzle 
diameters rather than the relaxed jet diameters were used. 
This linear dependence corresponds to the prediction for a 
low viscosity liquid and yet their apparent viscosities at 
low rates of shear are not indicative of such behavior. If, 
however, the MPA dispersions have a structural relaxation 
time which exceeds the lifetime of the jet then the apparent 
viscosity maintained m the jet will be approximately that 
existing at the capillary exit. Such a fluid would propagate 
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disturbance waves for an id'eal liquid .- 

The breakup length of capillary jets of 0.2 % Carbopol 
increases with -the nozzle diameter and with the average 
fluid velocity, as shown m figure 6-29 . At high veloci- 
ties, data for the 0.0216 cm. nozzle indicate that the 
breakup length becomes independent of the nozzle diameter. 

No maximum is observed m the curves within the experi- 
mental range of velocities investigated, as is seen for 
Newtonian and 0.1 % Carbopol jets. The breakup length 
data do not appear to pass through the origin and the 
reason for this is felt to be due to the unrelaxed molecular 
structure as explained previously for the MPA 60 dispersions 
(figure 6-27) . Similarly, the breakup lengths are not 
indicative of any exceptional stability which would be 
expected if the 0.2% Carbopol had regained its 'at rest' 
gel structure after profile relaxation. 

A crossplot of the figure 6-29 data, taken at a 
laminar velocity of 500 cm/sec is given m figure 6-30 and 
shows the dimensionless breakup length- L/2a 0 to vary 
linearly with the nozzle diameter. This dependence does not 
correspond to either of the Newtonian predictions for low 
and high viscosity liquids and this behavior is not yet 
understood . 

Solutions of 0.6% Carbopol are semi-rigid gells with 
very high apparent viscosities at low rates of shear, as may 
be seen from the data of figure 4-8 . The classical 
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stability theory of Weber, which assumes an initially 
relaxed liquid column at its zero shear viscosity, would 
predict complete stability for capillary jets of 0.6% 
Carbopol. Figure 6-31 shows that 0.6% Carbopol jets are 
stable at low velocities but are increasingly unstable at 
high velocities. This behavior is reminiscent of capillary 
jets of 12.4% silica dispersions (see figures 6-8 through 
6-10) . Again, m order to explain the data it must be 
assumed that the liquid column is not initially relaxed at 
its zero shear viscosity, but instead is substantially at 
the apparent viscosity possessed at the capillary exit. 

Thus, only at the very high velocities used to obtain the 
figure 6-31 data is the apparent viscosity low enough for 
the 0.6% Carbopol jet to be destroyed. 

The specific dependence exhibited m figure 6-31 
cannot be explained m terms of a shear dependent viscosity 
whose power law exponent is less than unity. However, at 
these very high velocities, air resistance is significant 
and it has been shown for Newtonian and 0.1% Carbopol jets 
under similar conditions that a decreased stability with 
increasing velocity is also displayed. It is likewise 
possible that, at these very high shear rates, some irrever- 
sible structural degradation might be occurring. 

In order to give some interpertation to the experi- 
mental data, the breakup lengths for capillary jets of 0.2% 
and 0.6% Carbopol solutions were compared to the Newtonian 
expectation that 
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Figure 6-31 BreaKup length of a capillary jet of 0.0}% Cartoopof as a 
function of fluid velocity; nozzle diameter- 0.0414 cm. 
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L/2a^ = C(We^ + 3 We^/Re^) (2.9) 

Two extreme cases were considered for the propagation 
of disturbance waves m viscoinelastic jets. The differences 
in the models will be reflected m the numerical value of 
the Reynolds number used m equation (2.9) . The first 
assumes that, outside the capillary, the fluid 'instantly' 
relaxes from all stresses generated during the capillary 
flow and also 'instantly' regains the gellified structure 
possessed under zero shear conditions. Since the shear 
rates associated with the growth of small, symmetrical 
disturbances are insufficient to cause significant departure 
from Newtonian behavior, the disturbance waves are assumed 
to propagate at a rate corresponding to the zero shear 
viscosity. However, the viscomelastic test fluids have 
a yield stress and do not display any clearly defined zero 
shear viscosity. For compantive purposes, the apparent 
viscosity at a shear rate of 0.1 sec"- 1 - will be used to 
approximate the zero shear viscosity. This may be obtained 
from figure 4-8 . 

In the other case considered, it was assumed that the 
time required to reform the .zero shear structure of the 
fluid is of the same order as or larger than the" lifetime 
of the jet. The liquid column, in this instance, is also 
taken to be initially relaxed from all prior stresses. 
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However „ because of' the long structural formation time, 
the viscoxnelastic fluid behaves as a Newtonian fluid with 
an average apparent viscosity approximately equivalent to 
that possessed at the nozzle exit* This behavior would be 
reflected in the use of a generalised non-Nevtonian 
Reynolds number such as that de fanned by Metzner and Park 
(53), Re^ = (2a 0 ) n (V 0 ) 2 ” n „ The choice of this 

particular Reynolds number is not unique. All quantities 
are evaluated at the capillary flow conditions. The 
constants K and n are fluid parameters xn the power law 
relation 


^^12, (1,14) 

Table 6-2 compares the breakup lengths predicted for 
each of these two proposed models with representative 
experimental data for capillary jets of 0,2 % Carbopol and 
0,6 % Carbopol o It is evident from the results that the 
assumption of a visco inelastic jet propagating disturbances 
at a rate corresponding to its zero shear viscosity is 
totally inadequate for these fluids. The predicted breakup 
lengths for this model are many orders ^of magnitude higher 
than the- experimental values. But the hypothesis that the 
viscosity possessed xn the capillary xs maintained through- 
out the lifetime of the jet provides results that are 
generally -consistent with those ,measured,A similar 
comparison with the MPA 60 dispersions was not considered 



172 


feasible because the thixotropic nature of these fluids 
precluded the determination of the apparent viscosity 
at low shear rates 0 Data obtained with capillary jets 
of 0.1 % Carbopol -was also excluded because of its unexplam 
able enhanced instability compared to invxscxd jets. 

Encouraged by the general agreement of the experimental 
data for the Carbopol solutions in table 6~2 with the 
assumption of a structural formation time whicn is large 
compared to the jet breakup time, breakup length measure” 
ments for all the Carbopol solutions and MPA 60 dispersions 
are plotted in the form of 


L 

2(X 0 




where Re 1 as given by 



(6.32) 


(6.32a) 


(The results are shown in figure 6<=32 where the solid lines 
indicate the range of the Newtonian breakup length data 
determined from this work and the measurmen ts of Grant 
(33) and Kroesser (39). General agreement is indicated 
between the vxscoinalastic and the Newtonian data when 
plotted in the form of equation (6.32) . This accord 
represents a strong statement for the belief that the 
gel structure is not formed td any significant extent 



Table 6-2. Comparison of calculated breakup lengths based on assumed 

viscosities with experimental values for viscomelastic nets. 


Fluid 

Nozzle 

Diameter 

Average Fluid 
Velocity-Vo 

(L/2ao) based 
on Yj , at 

(L/2a 0 ) 

expt 11. 

(L/2a 0 ) based on 
V at &=‘Q.l sec -1 

0-2% 

0.0216 cm 

300 cm/ sec 

.... (^i 

120 

120 

45,000 

Carbopol 

0.0263 

300 

111 

137 

45,000 


0.0414 

300 

133 

215 

45,000 


0.0868 

300 

182 

380 

45,000 

0.6% 

0.0414 

3700 

1610 

4700 

10 7 

Carbopol 


4900 

1880 

1000 

io 7 
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within the lifetime of the jet, The same conclusion was 
reached from an examination of the profile relaxation 
data presented in chapter 5® Had the zero shear viscosity- 
been used f instead of the apparent viscosity at the shear 
rate within the capillary, the data points would _have 
fallen several orders of magnitude away from the Newtonian 
data. 

The experimental data for the viscoinelastic jets 
should not be expected to correlate to a high degree for a 
variety of reasons. Undoubtedly, there is some viscosity 
buildup during the lifetime of the jet so that the use 
of the average apparent viscosity at the nozzle exit repre- 
sents an idealization. Secondly, the use of equation (6,32) 
was not derived from any fundamental consideration of the 
equations of motion but is merely an arbitrary modification 
of the Newtonian relation, equation (2 0 9), Considerations 
as to the effect of the radially varying, apparent viscosity 
actually existing at the nozzle exit have been omitted 
from the arguement. The constants K and n , determined 
from the capillary viscometer measurements, are not known 
with the same degree of accuracy as is a Newtonian viscosity. 

It is not uncommon that literature .values for these constants, 
determined by different investigators > differ by + 30 % , The 
magnitude of . the generalized non-Newtonian Reynolds number- 

t 

Re 0 is especially sensitive to the value of the exponent- n . 
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Finally, the possibility of shear degradation or of ultra- 
violet light deterioration between the time of property 
measurement and the determination of jet breakup lengths is 
unaccounted for. 

The wavelength of the disturbance wave was measured 
from photographs in the region just prior to breakup. An 
apparent viscosity at that point was calculated from the 
Newtonian relation, equation (6.26), using the nozzle 
diameter instead of the relaxed jet diameter. Table 6-3 
compares this calculated viscosity at breakup with the 
apparent viscosity based on the average shear rate within 
the capillary. Differences between the two values should 
be indicative of structural formation within the liquid. 

The results are consistent with the findings of the 
breakup length data just presented that the structural 
relaxation time is of the same order as or exceeds the 
lifetime of the jet. Thus, the apparent viscosity for 
0.1% Carbopol based on wavelength calculations displays 
the same mviscid nature as the apparent viscosity deter- 
mined from capillary viscometer measurements. Capillary 
jets of 0.2% Carbopol exhibit disturbance wavelengths which 

i 

denote a limited viscosity buildup from that possessed at the 
nozzle exit. The apparent viscbsities calculated from wave- 
length data for 24% MPA 60 m xylene are generally below the 
apparent viscosities computed at the average shear rate 
within the capillary. Wavelength measurements taken from 
photographs for this fluid were not consistent and since 



Table 6-3. Comparison of apparent viscosities calculated from wave length 
data at breakup with apparent viscosities at the capillary wall 
shear rate for viscoinelastic nets. 


Fluid 

Nozzle 

diameter 

Average Fluid 
Velocity-Vo 

Wave length 
at breakup 

Apparent 
Viscosity from 
eq. (6.26) 

Apparent 

Viscosity 

MPA60 

0.0216 cm 

328 cm/sec 

0.11 

cm 

0.10 poise 

0.35 poise 

in xylene 


474 

0.12 


0.15 

0.31 



630 

0.12 


0.15 

0.29 


0.0263 

382 

0.25 


1.7 

0.35 


0. 0.414 

286 

0.23 


0.21 

0.43 



448 

0.19 


0.03 

0.38 



562 

0.20 


0.07 

0.37 

0.2% 

0.0216 

228 

0.17 


0.79 

0.08 

Carbopol 

0.0263 

645 

0.20 


0.75 

0.058 


0.0414 

472 

0.30 


0.85 

0 . 08 

0.1% 

0.0414 

300 

.15 


0 

— 0.01-0.03 

Carbopol 

» - 

440 

.15 


0 

—*0.01-0.03 

0. 0868 

310 

.35 


0 

y~-* 0.01—0.03 



410 

.35 


0 

-'0.01-0.03 


177 
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changes of only 0.1 cm. significantly affect the calculated 
viscosities, the results are not considered accurate. 
However, all the viscomelastic test fluids showed terminal 
viscosity values that were of the same order as the appar- 
ent viscosities possessed at the capillary exit. The model 
of a viscoinelastic jet propagating disturbances at a 
rate corresponding to its zero shear viscosity is clearly 
invalid since the wavelengths demanded by this model are 
several orders of magnitude larger than the experimental 
values. 

6.6-3 CAPILLARY JETS OF VISCOELASTIC LIQUIDS 

In Newtonian jets, a disturbance is propagated as an 
exponentially growing wave with a constant wavelength. 

There is almost no difference between the place where the 
wave is strongly evident and, the place of final breakup. 
Capillary jets of 0.25% SCMC, a very slightly viscoelastic 
liquid, initially propagate disturbances m a similar 
manner but, as seen m figure 6-11, the growth of the wave 
is arrested and a string of droplets connected by threads, 
which thin with distance, is formed. Thread formation will 
be seen to be a characteristic of all viscoelastic jets 
used in this investigation and the reasons for its presence 
are not yet known. Final breakup occurs much later and, m 
this region of non-linear disturbances, the predictions of 
linear theory are inapplicable. 

The breakup length for 0.25% SCMC jets is therefore 
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defined as the distance from the nozzle where the disturb- 
ance wave amplitude is comparable to the original jet 
radius, a definition which would give satisfactory results 
for Newtonian jets, using this criterion, the breakup 
length and wavelength are found to lie between those 
predicted for an mviscid jet and those for a Newtonian 
fluid of the same zero shear viscosity, as shown m figures 
6-33 and 6-34 respectively. This result is m agreement 
with the theoretical predictions. The nozzle diameter, 
rather than the relaxed jet diameter, was used m the 
calculations as information on the latter is not available 
for capillary jets of 0.25% SCMC. 

Viscoelastic jets of 0.25% Separan and 0.25% Polyox, 
shown in figures 6-1 and 6-12 ^ possess a unique initial 
disturbance distance, which represents the distance from 
the nozzle to the sudden appearance of a large droplet. 

Values for this distance are not as reproducible as for the 
Newtonian jet breakup length and were found to vary by i 
30%. The initial disturbance distance for a capillary jet 
of 0.25% Separan increases with the fluid velocity and the 
nozzle diameter^ as shown in figure 6-35. A crossplot of 

I 

the data, taken at a fluid velocity of 400 cm/sec, which is 
within the linear range, is given m figure 6-36 and illus- 
trates the strong dependence of the (initial disturbance 
length/nozzle diameter) ratio on the nozzle diameter. As 
the nozzle diameter increases, the fluid becomes less sheared 
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and the (initial’ disturbance length/nozzle diameter) 
ratio becomes independent of diameter. This result would 
be expected for a highly shear sensitive liquid. 

The breakup length of capillary jets of dilute (0.05%) 
and moderately concentrated (0.25%) Separan, based on the 
point where the threads break, is shown m figure 6-37 as 
a function of the fluid velocity and the nozzle diameter. 
Determination of this distance was made by visual obser- 
vation of the liquid column under stroboscopic light and 
the accuracy of these measurements therefore dep&nds upon 
the ability to detect the finest threads. Values are found 
to be reproducible to + 10%. The breakup length of Q.05% 
Separan jets increases linearly with velocity and, at high 
velocities, becomes independent of it. Data points for all 
three nozzle diameters form a single continuous curve indi- 
cating that the breakup length- L/2a 0 is not a function of 
the capillary diameter. The breakup length of 0.25% 
Separan jets is seen to increase with fluid velocity and 
nozzle diameter. At sufficiently high fluid velocities, 
the breakup length also becomes independent of fluid velo- 
city. This behavior was exhibited for several of the vis- 
comelastic jets. 

The 'droplet-thread' configuration is the form 
eventually assumed by both the dilute and the moderately 
concentrated viscoelastic jets. Threads are seen to thin 
with increasing distance from the nozzle while the droplet 



184 



Figure G-37 BreaKu-p length of a ca.pUla.rv/ jet of 0.05% 
Separan and. 0.25% S eparan as a 
function of the average f Juid. velocity 
and. the nozzle diameter; 
nozzle diameter- 0-0.0'2I6 cm , 
A-0.02G3cm.j X- 0.0414 cm., □-0.0868cm., 
0,05% Separan, 0.25% S eparan. 
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diameter increases „ By photographing a 0.25% Separan jet 
along its length until breakup occurred, the droplet dia.- 
meter and thread length were determined as a function of 
distance along the liquid column. The computed results 
are based on 20-30 droplets and thread lengths at each 
location. Original data may be found in appendix B 

Figure 6-38 is a plot of the length mean droplet 
diameter vs. distance from the nozzle and shows that the 
mean diameter increases until about halfway through its 
flight and thereafter remains constant. Also shown m the 
same figure are the Sauter mean diameter and the weight 
mean diameter which are defined as follows 

D s = (6.33.) 

D w = ^D x 4 / SD i 3 (6.34) 

where D s and D w are the Sauter mean diameter and the weight 
mean diameter, respectively. The fact that all three 
measures of average diameter show the same trend indicates 
that the distribution remains similar. The standard devia— 
tion at each location is between 25-35%. 

The thread lengths between droplets were found to be 
widely distributed and showed no discernible pattern with 
increasing axial distance. At any one location, the values 
ranged from fractions of a centimeter to 6 cm. with an 
average of 2.5 cm. Thread lengths for 0.05% Separan jets 
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Figure 6*3 S Growth of droplets on a capillary jet 
of 0.Z5 /£> Separan as a function of the 
axial distance from the nozzle; 
nozzle diameter- 0.0216 cm., average 
fluid velocity- 1050 cm /sec., Q-lengf b 
mean diameter, A-5a.uter mean 
diameter, O- weight mean diameter 
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varied in a similar manner except that the threads are 
shorter. 

The change in droplet diameter with increasing axial 
distance and the random nature of the thread length at any 
one location both indicate that dynamic phenomena are occur- 
ring along the liquid column. Through the use of high 
speed motion pictures / taken at 3000-7000 pictures per 
second m the region where the droplet diameter is still 
increasing, the nature of these changes was investigated. 

Selected frames are shown m figure 6-39 for a 
capillary jet of 0.05% Separan and m figure 6-40, photo- 
graphs (a)-(f), for a 0.25% Separan jet. The droplets are 
seen to move at different velocities along the threads which 
create collisions m which new droplets of larger diameter 
are formed. Smaller particles decelerate more rapidly than 
larger particles due to frictional drag and this would lead 
to collisions of the type observed. However, this does not 
explain why the viscoelastic jet prefers the ' droplet- thread 1 
configuration m the first place. One significant observation 

i 

resulting from an analysis of the motion pictures is that 
the liquid column, far from being relaxed, continually 
undergoes large amplitude, axial deformations which require 
the existence of tensile stresses for their propagation. 
Photographs (g)-(j) of figure 6-40 show an instance where 
a whole section of the jet overtakes another section which 
was previously m front of it and then passes it. This 
phenomena is not as yet understood. 
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While thread thinning data for all the viscoelastic 
jets as a function of axial distance show considerable 
scatter, a dimensionless plot of the fraction of the 
original radius- a'/cio vs. the fraction of the breakup 
length- L'/L has the form shown m figure 6-41. Thus, 
breakup of the jet by disruption of the thinnest thread is 
imminent when the thread diameter is approximately % of the 
original diameter. Just before breakup occurs, the threads 
thin quite rapidly. In this region, the stability of the 
threads is quite remarkable. For example, m the portion 
of the liquid column immediately proceeding the breakup 
point, the threads connecting a series of droplets m a 
0 .,05% Separan jet (with 2aQ=0.0216 cm and ¥^=460 cm/sec) 

r 

decreases m thickness to 0.004 cm from 0.01 eta. over a 
length of 8 cm. A Newtonian liquid of the same zero shear 
viscosity would breakup in only 1-2 cm. under these 
conditions. 

The average diameter of the droplets formed after 
breakup as well as the size distribution was determined for 
capillary jets of glycerin, 0.05%~Separan and 0.25% Separan. 
Data for the Newtonian fluid was obtained m order to 
compare the accuracy of the measuring technique with 
published correlations for these liquids. Glycerin, m 
particular, was chosen because its Newtonian viscosity is 
approximately the same as the zero shear viscosity of 0.25% 
Separan. 




FRACTION OF ORIGINAL RADIUS- 
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Figure <3-41 Correlation of all viscoelastic test 
data for the decrease in jet radius 
with increasing axial distance; 

critical rad. taS at inception of 

breaKup. 
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Duffie and Marshall (17) investigated the droplet size 
distribution to be expected from the breakup of low speed, 
low viscosity Newtonian jets using capillary diameters and 
fluid velocities similar to those employed m this study . 
They found that the average, geometric mean diameter-Dg for 
all fluids investigated could be correlated by the empirical 
equation 

0.56 -0.10 

D g = 36 (2a Q ) (Re 0 ) (6.35) 

where Dg and a Q are m microns. D^. is defined as 

Ln (D g ) « ^Ln^) A (6.36) 

and N is the number of observations . There is some 
uncertainty whether equation (6.30), which was obtained with 
low viscosity liquids, can be extended to high viscosity 
fluids such as glycerin. 

The experimental data for glycerin droplets at 30°C 

i 

are compared to the predictions of equation (6.35) in 
table 6 -4 and the generally good agreement between the two 
indicates that the accuracy of the Pleasuring technique is 
comparable to those of other investigators. Data taken at 
a point 4 cm. past breakup shows a larger droplet diameter 
than the corresponding data taken at breakup and is due to 
droplet collisions and subsequent coalescence. A signi- 
ficant discrepancy does exist from the predicted value for 
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droplets measured at the lowest fluid velocity. This 
resulted from the jet's tendency to break into segments, 
rather than discrete droplets, with an accompanying 
increase m the particle size eventually formed . 


Table 6-4. Comparison of experimental droplet diameters 
for glycerin with tEe predicted values 'from 
equation (6.35). 

nozzle diameter-0.0216 cm. 
fluid temperature- 30°C. 


Experimental geometric 
mean droplet diameter 


Average fluid at 
velocity breakup 


4 cm. from Predicted values 

breakup - eg. (6.35) 


403 cm/sec 0.089 cm. 0.094 cm. 


0.070 cm. 


678 


0.069 0.075 


0.066 


937 0.071 0.072 


0.064 


The results for the average droplet diameter of the visco- 
elastic jets and of glycerin are summarized m table 6-5. 
Data taken for the glycerin and 0.05% Separan droplets 
showed a coefficient of dispersion (standard deviation/ 
length mean diameter) of 14-21% while that of the 0.25% 
Separan was 39%. The disparity arises from differences m 

I 

the nature of their breakup. From photographs taken at the 
point of breakup, it was observed that, m the dilute 
viscoelastic fluids, the fine connecting threads are 
absorbed into the droplets and a relatively uniform size 
distribution is produced. The comparatively small sepa- 
ration between the length mean diameter, the Sauter mean 



Table 6-5* Droplet diameters produced from the breakup of 
_ Newtonian and viscoelastic -jetst. 


Fluid 

Nozzle 
diame ter 

V o 

d l 

D 

s 

D 

w 

Coefficient of 
Dispersion 

0 . 05% 

0.0263 cm 

222 cm/sec 

0.106 cm. 

0.115 cm. 

0.119 cm. 

19.8 % 

Separan 


312 

0.094 

0.104 

0.111 

22.3 



435 

0.090 

0.100 

0.105 

24.4 



490 

0.083 

0.091 

0.095 

21.6 



715 

0.083 

0.088 

0.090 

16.8 



937 

0.084 

0.088 

0.090 

15.5 



1473 

0.080 

0.087 

0.091 

21.3 


0.0414 

349 

0.113 

0.121 

0.125 

17.8 



420 

0.125 

0.130 

0.133 

14.4 



505 

0.121- 

0.126 

0.129 

14.1 



722 

0.122 

0.130 

0.133 

17.2 


0.0868 

723 

0.2 22 - 

0.238 

0.245 

18.9 

0 . 25% 

0 . 0216 

1066 

0.133 

0,166 

0.174 

39.0 

Separan 





- 


Glycerin 

0.0216 

403* 

0,090 

0.096 

0,099 

17.8 



403** 

0.096 

0.103 

0.106 

18.8 



678* 

0.070 

0.076 

0,078 

20,0 



678** 

0.076 

0.082 

0.086 

21.0 



937* 

77 7C 

0,073 

0,078 

P.080 

20 6 



937 

0.073 

0.078 

0.081 

19.2 


Data taken at breakup 

Data taken 4 cm. past breakup 
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diameter and the weight mean diameter reflect the homo- 
geneity of the dispersion. However, the threads m a 
0.25% Separan jet collapse and form small satellite 
droplets which widens the droplet size distribution. 

This is also apparent from the larger standard deviation 
and from the increased differences between the three 
measures of droplet diameter. 

It is seen from the table that the droplet diameter 
resulting from the breakup of a capillary jet of 0.05% 
Separan, issuing from a 0.0263 cm. nozzle initially 
decreases with increasing fluid velocity and then becomes 
independent of it. A comparison with figure 6-37 shows 
that both the breakup length and the droplet diameter 
reach constant values over the same range of fluid 
velocity. A constant diameter region is also exhibited m 
the data taken with the 0.0414 cm. nozzle. In this con- 
stant diameter regime, figure 6-42 shows that the dimen- 
sionless droplet diameter- D L /2a 0 , where Dj, is the length 

- 0.2 

mean diameter, is proportional to (2a 0 ) 

Equation (6 .30) predicts droplet diameters for the 
Separan solutions which are substantially below the 
experimentally observed values. This is expected since 

i 

the high speed movies given in figures 6-39, and 6-40 have 
shown that the droplets formed on the intact capillary jet 
are continually colliding and coalescing to create larger 
diameter droplets. 
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Figure The effect of rcz.z.le diameter or 
the droplet diameter for a 
capillary jet of 0.05% Separan; 
average fluid velocity -17.2 cyr\/ssc. 
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The droplet size distributxon of glycerxn and of 0.d5% 
Separan follow a normal distrxbution as shown xn figures 
6-43 and 6-44, respectively, whereas the results of 0.25% 
Separan, given xn fxgure 6-45, devxate from thxs pattern. 

As explaxned previously, thxs occurs because the capillary 
jets of 0.25% Separan form small satellxte droplets at 
breakup along wxth the maxn droplets. Thus, there wxll be 
a preponderance of these small droplets xncluded and the 
sxze dxstrxbutxon wxll assume the shape of a bimodal dxs- 
trxbutxon. This xs more easxly seen xn the hxstograms of 
figure 6-46, which were constructed from representatxve 
droplet data for each of the three test fluxds. 

6.6-4 DISCUSSION OF EXPERIMENTAL RESULTS 

Fluxds whxch exhxbxt a yxeld stress or a shear depend- 
ent vxscosxty possess a three dimensional structure that 
requxres some fxnxte time to form or break down upon removal 
or applxcatxon of shearxng stresses. Thxs relaxatxon txme, 
xn some cases, could be longer than the txme for the 
velocity profxle and the shear stresses to reach equxlx- 
brxum. Thus, the hydrodynamxc behavxor of such a non- 
Newtonxan fluid should be related to at least two dxfferent 
relaxation txmes, a structural relaxatxon txme and a 
stress relaxatxon txme. Whether or not these relaxatxon 
txmes are sxgnxfxcant depends upon the txme scale of the 
experxment. A jet breakup experxment xs xdeally suxted to 
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Figure (o-Ar'b Droplet size cUsiribuf ion fora 
capillary jet of glycerin *, 
nozzle diameter - 0,0216 cm , 
O-data fcxKen 36 cm. from 
noziie, G-d.ata taKen 40 cm. 
from nozzle. 
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CUMULATIVE PERCENTAGE 


Figure <5-4-4 Droplet size cUstnfou/fion fora capillary 
jet of 0.05% 6 eparan as a f unction of the 
average f lu,?d. velocity and. the nozzle 
diameter* nozzle diavneter-A-0.02G3cm 3 
X-0.04-14 cm, D-0,08G8cm. 
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Figure <5-44 continued 
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Figure G-45 Droplet size distribution for a 
capillary jei of 0.25% S epa.ran- 
nozzie cliamefer-0.02IG cm, 
average f luict velocity-1045 cm/sec. 
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DROPLET D1AM£T£R-[Dl/'2oio] 

Figure 6-46 Comparison of droplet size distribution 
for cap! Mary jets of glycerin, 0.05% 
Separan and, 0,25%Separan - 

nozzle average 

fluid. cti ameter fluid, velocity 

glycerin 0.02i6cm R37 cnn/sec 

0.0 5% Separan 0,0263 S37 

0.25%Separan 0,0216 1045 
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demonstrate the effects of relaxation phenomena which are 
of the same order as the lifetime of the liquid column, 
about 0.001-0.01 seconds. 

The experimental results of this paper showed that the 
viscomelastic test fluids could all be atomized readily. 

This is surprising since several of them- the 0.6% 

Carbopol, the MPA 60 dispersions and the aqueous silica 
dispersion are highly gellified fluids when at 'rest. In 
particular, the silica dispersion, whichhas a relatively 
high yield stress, broke up m a manner similar to a low 
viscosity capillary jet m turbulent flow. With the 
exception of 0.6% Carbopol and the silica dispersion, the 
disturbance in the non-Newtonian inelastic fluids was 
propagated as a growing wave with a well defined wave- 
length. The two aforementioned fluids broke up only at 
very high velocities and, under these conditions, had an 
irregular, non- axisymme trie appearance. 

Classical Newtonian theory assumes the capillary jet 
to be initially relaxed from all prior stresses and at its 
zero shear viscosity. This model predicts complete stabi- 
lity for the viscomelastic test fluids since the small shear 
rates associated with the growth of infinitesimal, symmet- 
rical disturbances would not cause significant departure 
from Newtonian behavior. In the theory outlined earlier 
in the paper, a capillary jet of an initially relaxed 
viscoelastic liquid with a yield stress could be stable. 
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if the elastic modulus were such as to make the complex 

\ 

viscosity unbounded. 

For either of these theories to be applicable/ the 
structure which gives rise to the yield stress or the 
shear dependent viscosity would have to be formed within 
a time period which was short compared to the jet breakup 
time. It has been shown in the literature, e.g. Carver 
and Van Wazer (7) , that certain fluids require several 
minutes or hours to completely regain- their 'at rest* 
structure. Fischer et.al. (21) have presented data that 
demonstrates how the level of shear, to which the material 
was previously subjected, also influences the structural 
formation time . This latter point might explain why the 
silica jet was stable at the lower velocities or shear 
rates and atomized readily at the higher velocities. 

The breakup length and wavelength data for the visco- 
inelastic jets cannot be explained in terms of a iic^uid. 
column that relaxes instantly to regain its structure and 
then propagates disturbances at a rate corresponding to 
its zero shear viscosity. Such a model predicts results 
that are several orders of magnitude higher than experi- 
mentally observed values. 

The measurements are reasonably consistent, however, 
if it is assumed that all prior stresses relax shortly 
beyond the capillary exit but the time required to reform 
the gellified structure is of the same order or exceeds the 
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lifetime of the jet. In this model, the liquid column 
behaves as a Newtonian liquid whose viscosity is approxi- 
mately the same as the apparent viscosity of the visco- 
melastic fluid, evaluated at the average shear rate within 
the capillary. The proposed behavior is also m agreement 
with the experimental observation that the breakup length 
vs. fluid velocity curves do not appear to pass through 
the origin. However, no explanation is available for the 
increased instability of capillary jets of 0.1% Carbopol 
over an mviscid fluid, when both are m laminar flow. 

The linearized stability analysis has shown that the 
growth rate of axisymmetric disturbances for any visco- 
elastic fluid, whose complex viscosity is bounded, is higher 
than that of a Newtonian fluid of the same zero shear 
viscosity. If the complex viscosity could readily be 
measured, then one could predict the growth rate aS a 
function of the wavelength, assuming that the liquid column 
is completely relaxed initially. However, for the dilute 
solutions used m these experiments, the complex viscosity 
cannot be determined and only a qualitative prediction can 
be made. 

In the experiments, growth rates were not measured 
directly and comparisons of jet stability were made on the 
basis of the breakup length. This approach gives good 
agreement with Weber's theoretical predictions for 
Newtonian liquids. In the region of linear wave propa- 
gation, the breakup length and the wavelength of a capillary 
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jet of 0.25% SCMC were found to be shorter than a Newtonian 
liquid of the same zero shear viscosity but greater than 
for an mviscid liquid. This result is m agreement with 
the predictions of the linear theory. 

Contrary to what occurs m Newtonian liquids, the 
growth of the disturbance wave in 0.25% SCMC solutions 
does not lead to immediate breakup and the droplets remain 
connected for a considerable distance by thin threads. This 
thread formation seems to be typical of elastic liquids. 

No wave formation has been observed m highly elastic fluids 
such as 0.25% Separan and 0.25% Polyox. Even m rela- 
tively dilute solutions (0.05% Separan), the waves which do 
appear initially are damped out. Instead, disturbances 
propagate as isolated droplets connected by random lengths 
of threads. In many cases, the threads are longer than any 
reasonable wavelength that could be associated with the 
droplet. These non-linear disturbances occur at a distance 
which is shorter than the breakup length of a Newtonian 
liquid of comparable zero shear viscosity. Thus, one cannot 
conclude that the viscoelastic jet is more stable with 
respect to the growth of axisymmetric waves. 

It has been observed for all the viscoelastic jets 
that the droplets which are formed are connected by threads 
which continually thin and lead to the eventual breakup of 
the jet. Two questions arise: why do the threads form and 
what causes their remarkable stability? While the original 



210 


jet does not show any increased stability as compared to a 
Newtonian jet, the connecting threads have a lifetime 
several fold higher than expected. Viscoelastic jets 
extruded from a nozzle which is of the same order of thick- 
ness as the threads do not exhibit this increased stability. 

Two reasons might explain this difference. First, a 
viscoelastic jet extruded from a nozzle which is of similar 
thickness as the threads has just experienced a high 
shear rate. If the breakup time and the structural 
relaxation time are of comparable magnitude, the jet will 
break up before the fluid has regained its zero shear 
viscosity. However, the thin thread at a distance from the 
nozzle is derived from a much thicker jet, which was less 
severe^ sheared in the capillary and which has also had time 
to partially regain its zero shear viscosity. The second 
and more important difference is that the thread is under 
some stress, as shown by the fact that it continues to thin. 
This stress might be caused by either the surface tension 
forces at the connection to a droplet or by the relative 
motion of two droplets. Looking at photograph (a) of figure 
6-14, one might question this explanation as the liquid 
column looks quite relaxed. However, the high speed movies 
indicate that the threads and droplets are constantly chang- 
ing their configuration and are undergoing large scale 
deformations, as seen m figures 6-39 and 6-40. 

The stability of a viscoelastic thread has similarities 
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to the unique ability of viscoelastic liquids to be spun 
into fibers. During thread thinning, the fluid undergoes 
a flow which is characterized as a stretching or extrusion 
flow. It is well known that the behavior of elastic 
liquids under tension is completely different from that in 
shear flow m the sense that viscoelastic liquids possess 
an elongational viscosity which increases with the 
stretching rate. Thus, elastic liquids of low viscosity 
can be spun into fibers whereas Newtonian liquids are 
SpinnabLe only at high viscosities. 

i 

6.6-5 CONCLUDING REMARKS 

The results presented m this chapter have shown the 
limitations inherent in using assumptions, which are valid 
for capillary jets of Newtonian liquids, m the analysis 
of non-Newtonian fluids without prior experimental 
verification. In particular, it is evident that stress 
relaxation phenomena occuring outside the capillary exit 
can play a significant role m the ultimate behavior of the 
liquid column. The atomization of a non-Newtonian jet thus 
involves considerations of the fluid's transient and steady 
state flow properties. 

A very practical result of this research is the finding 
that highly gellified liquids are unstable if their 
structural relaxation time is of the same order as or 
exceeds the lifetime of the jet, about 0.001-0.01 seconds. 
Additives such as Carbopol, m a 0.1% concentration, have 
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the ability to thicken a solution yet are able to enhance 
the atomization. The classical notion that viscosity always 
increases the stability of a jet must be re-examined in 
light of these data. 
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CHAPTER 7 

SUMMARY AND CONCLUSIONS 

The objective of this research has been to carry out 
an< experimental program which would provide an initial 
understanding into the behavior of capillary jets of 
non-Newtonian liquids. 

Two categories of fluids were chosen for study; visco- 
inelastic materials and viscoelastic liquids. The viscomr 
elastic fluids were comprised of a wide variety of substances 
and included organic polymer solutions, an inorganic 
dispersion and organic polymer dispersions. They ranged m 
consistency from relatively thin liquids to stiff gels and 
could be characterized as possessing a shear dependent 
viscosity, a yield stress and an absence of elastic phenom- 
ena. The 'Viscoelastic solutions were all composed of 
organic polymer solutions of a class which exhibits marked 
drag reduction in turbulent flow. 

The experiments concentrated on profile decay measure- 
ments and on jet breakup studies. Data were obtained 
through high speed still and motion picture photographs 
taken along the length of the jet. 

Profile decay measumnaits refer to the change m jet 
diameter with axial distance as the emerging jet relaxes 
from the stresses generated during the previous shear flow. 

A linear perturbation analysis taken about the final relaxed 
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state showed that a Maxwell model fluid would relax faster 
than a Newtonian jet of the same zero shear viscosity. 

Experimental profile decay data obtained for Newtonian 
liquids confirmed the theoretical prediction that , within 
the range of variables investigated, the damping coefficient 
is a function only of the Reynolds number. Values of the 
damping coefficient were in fair agreement with the analytical 
solution. 

Profile decay data taken for a typical viscomelastic 
fluid- 0.6 % Carbopol could be satisfactorily correlated 
only by assuming that the structural formation time was 
large compared to the time required for stress relaxation. 

This model assumed that no viscosity buildup occured above 
that possessed at the capillary exit. This proposed explanat- 
ion was found to be consistent with wavelength and 
breakup length data for the viscomelastic jets For 
most of the test liquids, the breakup lengths corresponded 
to those for a Newtonian fluid of the same apparent 
viscosity as that possessed within the capillary. One 
notable exception was 0.1 % Carbopol which displayed 
an enhanced instability compared to water, when both were 
m the laminar range. Wavelength data taken near the point 
of breakup also confirmed the lack of any significant 
structural formation. All these observations lend support 
to the theory that the behavior of a capillary jet of a 
viscoinelastic liquid can be explained m terms of a 
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structural formation time which is large compared to the 
lifetime of the jet* 

A theoretical analysis due to Yerushalmx (94) had 
predicted that a fluid with a sufficiently high yield stress 
could lead to a completely stable jet. However, the develop- 
ment assumed that the net was initially relaxed at its zero 
shear structure. Such a conclusion would be inapplicable 
to the viscomelastxc liquids of this investigation which 
never reform their 'at cest ' structure within the breakup 
time of the 3 etc. 

A practical application of these experimental results 
is that highly gellifxed viscoxnelas tic fluids can he 
readily atomized . The stiffest liquid studied, 12 .4 % 
silica in water, had the appearance of a low viscosity 
jet m turbulent flow. It also appears that certain 
additives such as 0.1 % Carbopol are able to thicken a 
solution and yet enhance its atomization. 

Profile decay data for a typical viscoelastic fluid- 
0,,7 % SCMC showed decay rates slower than for a Newtonian 
jet of the same zero shear viscosity but faster than that 
for one of the same apparent viscosity,. This is in direct 
contradiction to the theory of profile relaxation 
discussed earlier. The results can be xnterperted as 
showing that some structural formation had occurred. Another 
possibility which can not be resolved by these data alone 
is that a slowly decaying axial tensile stress exists along 
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•the jet length which serves to retard the diameter change 
on a contracting column. However, high speed movies taken 
for several of the viscoelastic test fluids showed that 
large scale axial deformations are continually occuring 
and these require the existence of axial tensile stresses . 
The viscoelastic jets were seen to propagate disturbances 
by means of a droplet-thread configuration and the 
preference for this unusual configuration is not yet 
understood „ Whi^e the evidence can not be regarded as 
conclusive, the data do point towards a limited viscosity 
or structural buildup within the lifetime of the jet and 
do demonstrate vividly the effects of axial normal stresses. 

The linearized stability analysis for the breakup of 
a viscoelastic jet had predicted a higher growth rate of 
disturbances than for a Newtonian liquid of the same zero 
shear viscosity. This was partially confirmed for jets of 
0.25 % SCMC. But even in the most dilute fluids studied, the 
initial wave propagation of disturbances was arrested and 
evolved into the unique droplet-thread configuration of 
these fluids. 

While these studies have proposed some tentative 
theories to explain the behavior of capillary jets of visco— 
inelastic and viscoelastic liquids, they have also raised 
many questions and indicated the need for continued research 
in this area. 
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